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ABSTRACT 


Inpile  microwave  measurements  of  electron  density  in  neon-argon  and 
argon-cesium  plasmas  generated  by  fission  fragments  are  compared  with  values 
of  electron  density  predicted  from  a  reaction  kinetics  theory.  The  main  pur¬ 
pose  of  the  comparison  is  to  assess  the  validity  of  a  theory  designed  to 
describe  the  dominant  production  and  loss  processes  in  noble  gas-cesium 
thermionic  converters  over  a  wide  range  of  conditions .  For  the  neon-argon 
system  the  measured  and  predicted  values  were  in  good  agreement;  the  highest 
electron  density  was  1.0xl012cm"3  at  full  reactor  power  (<|>=1 .44xl013cm'2sec_1) . 
For  the  argon-cesium  system  the  agreement  between  theory  and  experiment  was 
less  satisfactory.  The  highest  measured  electron  density  at  full  reactor 
power  was  2.3xl012cm"3  —approximately  twice  the  computed  value,  and  this 
fflavlTnr.nl  in  the  curve  of  electron  density  vs.  the  ratio  of  Cs/Ar  occurred  at  a 
ratio  U>2xl0"^)  much  lower  than  that  predicted  (-o 3x10  ^) .  Furthermore  the 
electron  density  was  found  to  be  extremely  dependent  upon  the  temperature  of 
the  cavity  walls .  No  satisfactory  explanation  has  yet  been  found  for  this 
behavior . 

Favorable  electron  transport  properties  are  expected  to  make  the  fission- 
fragment  -generated  argon-cesium  plasma  a  good  candidate  for  use  in  a  nuclear 
thermionic  converter.  Theoretical  transport  properties  are  reported  for  this 
plasma  when  the  major  ion  loss  mechanism  is  ambipolar  diffusion  to  the  therm¬ 
ionic  diode  electrodes .  The  theoretical  model  is  being  expanded  to  include 
volume  loss  of  the  ions. 
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INVESTIGATIONS  ON  THE  DIRECT  CONVERSION  OF 
NUCLEAR  FISSION  ENERGY  TO  ELECTRICAL  ENERGY  IN  A  PLASMA  DIODE 

INTRODUCTION 


The  analysis  of  reaction  kinetics  and  electron  transport  in  a  noble-gas 
plasma  generated  by  fission  fragment  ionization  has  been  continued.  These 
studies  evaluate  such  a  plasma  for  application  to  a  nuclear  thermionic  energy 

converter.  (]_) 

The  rate  at  which  fission  fragments  produce  ions  in  noble  gases 

been  used  in  reaction  kinetics  studies  to  predict  the  steady-statylectron 
density  produced  in  both  a  neon-argon  and  an  argon-cesium  plasma.  In 
the  previous  reporting  period,' (2>  the  transport  of  electrons  across  a  fission- 
fragment -generated  plasma  was  considered  for  a  plasma  controlled  by  diffusion 
loss  of  the  ions,  and  the  description  of  an  inpile  device  was  presented  for 
the  measurement  of  the  electron  transport  properties  of  such  a  plasma. 

Tn  this  report  we  present  first  the  experimental  results  of  inpile  micro¬ 


wave  experiments  to  measure  the  electron  density  produced  in  both  neon-argon 
and  argon-cesium  plasmas.  These  results  are  compared  with  the  values  pre¬ 
dicted  by  the  reaction  kinetics  theory.  In  the  second  section  we  present 
the  results  of  parametric  studies  on  electron  transport  in  a  diffusion  con¬ 
trolled  plasma.  An  initial  transport  device  was  operated  inpile  but  the 
results  to  date  are  preliminary. 


OBJECTIVES 


The  objectives  initially  set  for  the  current  reporting  period  were: 

1.  Continue  microwave  measurements  inpile  to  determine  experimentally  the 
electron  density  generated  by  fission  fragments  in  the  Ne-Ar  and  Ar-Cs 

plasmas . 

2.  Compare  measurements  of  electron  density  with  predictions  of  the  reac- 
tion  kinetics  analysis  under  various  conditions  . 
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3.  Operate  inpile  thermionic  diodes  incorporating  the  Ba0-U02-W  emitter  with 
argon-cesium  and  neon-argon  fillings  to  study  the  electron  transport 
through  the  fission-fragment-generated  plasma. 

4.  Extend  theoretical  work  on  electron  transport  through  noble-gas  plasmas 
to  include  volume  loss  processes  which  are  both  linear  and  quadratic  in 

electron  density. 


CONCLUSIONS 


1.  The  values  of  electron  density  in  the  neon-argon  plasma  (90  torr  neon, 
10“4  Ne/Ar)  as  measured  with  the  resonant  microwave  cavity  were  in  good 
agreement  with  the  values  predicted  by  the  reaction  kinetics  theory. 

2 .  The  highest  measured  electron  density  in  the  neon-argon  plasma  was 
1.0xl0^2cm”3  at  full  reactor  power  (2  MW)  and  (|)=1 .44x10  cm  sec 

3.  Variation  in  the  intensity  of  the  microwave  signal  to  the  cavity  (Ne-Ar) 
produced  little  effect  on  the  resonant  frequency  shift  since  the  cavity 
was  greatly  undercoupled  to  the  waveguide  at  the  higher  electron  densi¬ 
ties  where  some  effect  might  be  expected  • 

4.  Variation  in  the  cavity  wall  temperatures  (and  therefore  gas  tempera¬ 
tures)  also  had  little  effect  on  the  electron  density  for  Ne-Ar.  A 
much  larger  variation  was  predicted  by  the  reaction  kinetics  code. 
However  at  the  present  time  the  computations  include  the  temperature 
dependence  of  only  two  (diffusion  and  collisional  radiative  recombina¬ 
tion)  of  the  many  processes  involved .  The  temperature  dependence  of 

the  other  reaction  rates  are  being  studied  . 

5.  For  the  argon-cesium  microwave  cavity  (90  torr  argon),  a  variation  of 
Cs/Ar  from  10“^  to  10”3  was  obtained  simply  by  the  use  of  two 
nitrogen  circulating  tubes  to  remove  the  excess  fission  heat. 

6.  The  highest  electron  density  measured  for  Ar-Cs  was  2-3x10  cm  at  full 

reactor  power  (2MW)  (<|=1  Mxl013cm-2sec_1)  with  Cs/Ar  2x10  . 

7 .  The  reaction  kinetics  theory  had  predicted  higher  maximum  electron  den¬ 
sities  for  Ar-Cs  than  for  Ne-Ar  and  this  was  borne  out  experimentally. 
Even  so  the  maximum  electron  density  in  Ar-Cs  was  higher  than  that  cal¬ 
culated,  and  more  importantly,  it  occurred  at  a  much  lower  value  of 
Cs/Ar  (~>2xl0"5)  than  that  (~.  3x10“ 4 )  predicted. 
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8.  The  measured  electron  density  in  Ar-Cs  was  extremely  dependent  upon  the 
temperatures  of  the  cavity  walls .  In  one  run  an  increase  of  only  200  C 
(  0.017  eV)  in  the  average  wall  temperature  produced  an  order  of  magni¬ 
tude  increase  in  the  electron  density  for  the  same  corrected  value  of 
Cs/Ar  ratio. 

9 .  Explanations  of  this  anomalous  dependence  of  electron  density  upon 
temperature  have  been  sought  in  possible  surface  and  volume  reactions 
not  presently  included  in  the  reaction  kinetics;  so  far  none  appear 
capable  of  yielding  the  observed  temperature  behavior . 

10.  The  thermionic  electron  transport  model  with  a  diffusion-controlled 

12  -3 

plasma  is  applicable  to  electron  densities  below  2x10  cm  .  Above  this 
density,  the  volume  ion  losses  become  appreciable. 

11.  This  diffusion  model  with  a  uniform  volume  ionization  rate  from  fission 
fragments  predicts  considerable  variation  in  the  electron  density  between 
the  emitter  and  collector;  typically,  ne(max)/ne(min)~*40. 

12.  The  calculated  diode  current  was  very  sensitive  to  the  ion  generation 
rate  and  diode  spacing  as  expected.  For  the  plasma  limited  conditions 
studied,  the  diode  current  was  relatively  insensitive  to  the  emitter 
Richardson  current . 


FUTURE  PLANS 

1.  The  experimental  and  theoretical  study  of  reaction  kinetics  in  argon- 
cesium  plasmas  will  be  continued  in  order  to  establish  the  nature  of 
the  observed  critical  dependence  of  electron  density  upon  ambient  gas 
and  wall  temperatures . 

2.  Electron  transport  tubes  with  Ar-Cs  fillings  will  be  operated  to  measure 
electron  current  transport  as  a  function  of  electron  density  in  the 
plasma  as  determined  from  the  reaction  kinetics  study. 

3.  Theoretical  work  on  electron  transport  will  continue  with  the  aim  of 
developing  a  theory  which  will  satisfactorily  predict  current -voltage 
characteristics  and  performance  of  argon-cesium  plasma  converters. 

The  theory  will  be  tested  against  the  experimental  results  obtained  with 
the  electron  transport  tubes  . 
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4  .  Consideration  will  be  given  to  operating  a  transport  tube  with  an 

improved  collector  surface  in  order  to  increase  power  output .  Work  on 
the  influence  of  fission  fragments  on  these  surfaces  will  be  initiated. 
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SECTION  A 


INPILE  MICROWAVE  MEASUREMENTS  OF  ELECTRON  DENSITY 
IN  NEON -ARGON  AND  ARGON -CESIUM  PLASMAS 

ABSTRACT 

Inpile  microwave  measurements  of  the  electron  density  in  neon-argon 
plasmas  have  been  completed  and  are  reported  together  with  similar  measure¬ 
ments  for  argon-cesium  plasmas .  The  plasmas  were  generated  by  fission 
fragment  ionization  of  the  corresponding  mixed  gases  in  a  resonant  micro- 
wave  cavity.  Experimental  techniques  are  described  and  the  experimental 
data  are  analyzed  in  terms  of  the  reaction  kinetics  theory  which  predicts 
the  electron  density  from  the  various  computed  ion  source  and  loss  rates 
within  the  plasma.  For  the  neon-argon  plasma,  theory  and  experiment  agreed. 
With  argon-cesium,  higher  electron  densities  were  found,  as  predicted  by 
theory,  but  the  predicted  maximum  in  electron  density  occurred  at  a  much 
lower  concentration  of  cesium  than  had  been  calculated .  Also  a  temperature 
effect  was  observed  in  argon-cesium  whereby  the  electron  density  was  found 
to  increase  very  rapidly  as  the  temperature  of  the  walls  of  the  cavity  was 
increased  from  about  400  to  700°K  at  constant  cesium  concentration. 
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OBJECT 


The  objectives  of  this  study  were  (l)  to  measure  the  electron  density 
in  a  fission-fragment-generated  plasma  using  a  resonant  microwave  cavity 
and  (2)  to  compare  the  measured  electron  density  with  that  predicted  by 
the  reaction  kinetics  theory. 


CONCLUSIONS 


1.  The  measured  electron  densities  in  the  neon-argon  plasma  agreed  well 
with  theory. 

2.  The  variation  in  electron  density  with  gas  temperature  was  not  as 
great  as  expected  for  neon-argon.  However  the  temperature  depend¬ 
ence  of  only  the  diffusion  and  collisional  radiative  recombination 
rates  have  so  far  been  included  in  the  theory . 

3.  Use  of  the  new  wide  sweep  K-band  microwave  generator  permitted  much 

more  accurate  measurement  of  the  electron  density  and  disclosed  four 
more  resonant  peaks  (other  than  moc*e)  which  were  also  used  to 

measure  the  electron  density. 

4.  Microwave  measurements  on  the  argon-cesium  plasma  gave  higher  elec¬ 
tron  densities  than  for  neon-argon, as  had  been  predicted  by  theory. 

5 .  When  the  cesium  concentration  in  the  argon  varied  by  changing  the  cesium 
bath  temperature,  the  electron  density  passed  through  a  maximum  value 
but  at  a  cesium  concentration  much  lower  than  that  predicted. 

6.  Another  unexpected  (and  very  probably  related)  phenomenon  was  that 
the  electron  density  was  very  sensitive  to  the  temperature  of  the  cav¬ 
ity  walls  ( 400-700  °K) . 

7.  The  logarithm  of  electron  density  varies  linearly  with  the  inverse  of 
the  average  cavity  temperature  which  yields  an  "activation  energy"  of 
0.22  eV  for  an  unspecified  process. 

8.  Several  processes  additional  to  those  included  in  the  reaction  kinetics 

equations  have  been  briefly  considered  in  an  attempt  to  elucidate  the 
anomalous  temperature  variation  of  electron  density •  These  have  includ¬ 
ed:  (i)  the  time  for  the  cesium  to  diffuse  from  the  reservoir  into 

•  the  cavity  and  reach  an  equilibrium  value;  (ii)  a  delayed  equilibrium 
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of  cesium  arising  from  surface  effects;  (iii)  Langmuir  surface  ioniza¬ 
tion  of  cesium  on  the  cavity  walls;  (iv)  impurities;  (v)  interaction 
of  cesium-covered  surfaces  with  photons  from  the  plasma;  (vi)  volume 
ionization  of  cesium  by  argon  excited  states  which  are  not  metastable; 

(vii)  heteronuclear  ions;  and  (viii)  the  influence  of  energetic  elec¬ 
trons  from  the  Penning  ionization  process.  None  of  these  processes  act¬ 
ing  independently  in  the  present  experiment  appear  likely  to  be  the  cause 
of  the  observed  sensitive  dependence  of  electron  density  upon  the  cavity 
wall  temperature . 

9.  A  46  hour  run  with  the  Ar-Cs  cavity  at  high  temperature  and  high  electron 
density  showed  no  degradation  of  electron  density  with  time . 

INTRODUCTION 

The  noble  gas  thermionic  converter  with  the  plasma  generated  by  fission- 
fragment  ionisation  of  the  gas  has  been  described  1  and  theoretical  and 
experimental  studies  on  the  ion  generation  rate  in  such  plasmas  have  been 
reported .  ^  The  reaction  kinetics  analysis  of  mixed-gas  plasmas  and  pre¬ 
liminary  results  of  electron  density  measurements  in  Ne-Ar  were  presented  at 
the  San  Diego  meeting ;(3)  also  a  brief  report  of  similar  measurements  for 

(4) 

Ar-Cs  has  appeared  recently. 

In  our  theoretical  studies  of  the  reaction  kinetics  of  Penning-type  bin¬ 
ary  gas  systems  we  showed  that  considerable  enhancement  of  the  electron 
number  density  over  that  for  pure  gases  could  be  expected  by  the  proper  selec¬ 
tion  of  the  conditions  that  (l)  maximize  the  ionization  of  the  trace  gas 
atoms  by  metastable  excited  atoms  of  the  major  gas  species  (Penning  reaction) 
and  at  the  same  time  (2)  minimize  the  loss  rate  of  the  trace  gas  ion.  Reac¬ 
tion  rate  equations  for  the  atomic  and  molecular  ions  of  both  species  and  the 
metastable  atoms  of  the  major  species  in  the  plasma  were  solved  for  the  steady 
state  with  a  computer  to  find  the  influence  of  the  gas  pressure  and  composi¬ 
tion  on  the  electron  density.  The  results  of  these  computations  are  repro¬ 
duced  here  in  Fig.  1  for  neon-argon  and  for  the  argon-cesium  system.  The  gas 
temperature  (l300°K)  and  neutron  flux  (lxl013cm"  sec’  )  are  typical  conditions 
for  our  thermionic  converter  operating  in  a  research  reactor  (2  mW) .  In  both 
systems  a  singular  electron  number  density  maximum  was  found  at  a  pressure  of  . 
about  90  torr  and  at  a  trace  gas  concentration  of  about  0.01%.  It  was  to  verify 
this  predicted  behavior  that  the  inpile  microwave  measurements  of  electron  den¬ 
sity  were  performed.  ,  2 
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Fig.  1.  Contour  map  showing  the  singular  electron  number -density  maximum 
in  a  1300 °K  fission-fragment-generated  plasma  of  neon-argon  (a) 

-  ■  ‘  rth  of  the  container  is 
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and  argon-cesium  (b) .  The  diffusion  length  c 
1.6  mm  and  the  neutron  flux  is  lxlCr^cnT  sec- 
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EXPERIMENTAL  METHOD 


The  method  of  measuring  the  electron  density  (nj  in  a  plasma  from  the 
shift  ( AW)  in  the  resonant  angular  frequency  of  a  microwave  cavity  is  well 
known.  ^  Au»  is  related  to  ng  by 

(1) 

€0m  7  E2dV 

VJ 

where  E  is  the  electric  field  within  the  volume  V  of  the  cavity, is  the  fre¬ 
quency  of  the  empty  cavity,  e  and  m  are  the  electronic  charge  and  mass  of  the 
electrons,  €  is  the  permittivity  of  free  space  and  v  is  the  collision  fre¬ 
quency  of  the  electrons .  The  fission-fragment  ion  generation  rate  is  fairly 
uniform^  and  we  expect  to  operate  in  the  regime  where  the  volume  losses  pre¬ 
dominate  .  Thus  the  electron  density  is  assumed  to  be  uniform  over  V  and  the 
integrals  cancel  in  Eq.l  (after  removal  of  ng).  We  will  operate  at  pressures 
^100  torr  where  the  collision  frequency  is  of  order  v*^5xl0  sec  •  ■^ie 
resonant  angular  frequency  will  be  ~  1.5xl0118ec"1  (K-band  microwaves)  so  the 
factor  in  brackets  in  Eq.l  becomes  unity.  With  these  simplifications,  Eq.l 
reduces  to 

f 

n  (cm-3)  =  - -n  Af(GHz) 

e  4  X  10' 

where  f  is  the  resonant  frequency  of  the  empty  cavity  in  GHz. 
o 

In  order  that  the  microwave  cavity  could  be  mounted  inside  a  tube 
(i  .d  .=2-3/4  in.)  which  was  to  fit  a  standard  fuel  element  position  in^the  core 
of  a  swimming  pool-type  research  reactor,  we  selected  the  M020  mode  for  a 
right  circular  cylinder  for  the  cavity  geometry  as  shown  in  Fig.  2.  The  gen¬ 
eral  method  of  cavity  construction  was  similar  to  that  used  earliet  but  will 
be  repeated  both  for  completeness  and  since  Fig.  2  incorporates  improvements 
over  the  earlier  design.^ 

The  body  of  the  cavity  was  made  of  Kovar.  The  window  was  made  of  high  __ 
density  alumina  and  after  coating  with  molybdenum-manganese  it  was  brazed  to 
the  Kovar  retainer  with  copper.  A  short  section  (2-1/4  in.)  of  stainless 
steel  waveguide  was  brazed  to  the  Kovar  at  the  same  time  that  the  ceramic 
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TNV  2  Microwave  cavity  for  Ar-Cs  inpile  studies.  The  Ne-Ar  cavity 

was  identical  except  for  a  small  pump-out  tubulation  ins  ea 
rvf  +.Vip  cesium  bath. 


window  and  retainer  were  brazed  in  place .  The  inside  of  the  cavity  (less  the 
Kovar  cap)  was  electroplated  with  copper  (-0.004  in.).  Two  tantalum  discs 
(0.010  in.)  were  brazed  to  the  Kovar  cap  with  copper  and  the  U-235  foil 
(0.001  in.)  was  brazed  to  the  top  tantalum  disc  with  nickel.  After  the  Kovar 
cap  was  spot-welded  in  place  the  periphery  was  heliarc-welded  to  complete  the 
tube  construction.  This  cap  design  afforded  better  electrical  contact  and  gave 

better  signal  characteristics  than  the  earlier  design. 

me  M  mode  for  a  right  circular  cylinder  cavity  (Pig.  2)_was  selected 

with  f  =23.1  GHz  (K-band)  so  that  an  electron  density  of- 5x10  cm  would 
yield  a  frequency  shift &f-l  GHz.  This  mode  fixed  the  cavity  inside  diameter 
at  d=2 .272  cm  but  the  height  was  still  arbitrary.  For  the  two  cavities  dis¬ 
cussed  here  the  height  was  h=0.7  cm  (in  the  first  inpile  run  h=0.5  cm  ) . 

In  the  construction  of  the  first  He-Ar  microwave  cavity  considerable 
trouble  was  experienced  in  obtaining  a  vacuum-tight  seal  between  the  ceramic 
microwave  window  and  the  cavity.(1)  This  problem  was  solved  by  first 
bonding  the  ceramic  to  a  thin  walled  Kovar  tube  and  then  brazing  the  Kovar 
tube  to  the  cavity  wall  (Fig.  2).  Also  the  circular  ceramic  window  for  the 
first  cavity  was  made  as  large  in  diameter  as  possible  consistent  with  K-ban 
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waveguide  but  the  empty  cavity  was  undercoupled  nevertheless  and  severe  loss 
of  signal  resolution  developed  as  the  plasma  density  increased  (plasma  decreases 

f  cr  \ 

the  coupling) .  To  solve  this  problem,  for  the  cavities  reported  here,  the 
ceramic  window  was  changed  from  a  circular  to  oval  shape  in  order  to  give  a 
larger  cross  sectional  area .  This  change  gave  an  overcoupled  empty  cavity  and 
better  signal  resolution  with  a  plasma. 

After  construction,  the  cavity  was  baked  at  600°C  on  an  ultra-high  vacuum 
system*  and  evacuated  to  a  pressure  of  6x10  ^  torr.  Then,  under  vacuum,  the 
five  cesium  chromate  pellets  contained  in  a  glass  side-arm  were  carefully 
flashed  inductively  and  the  cesium  subsequently  flash  distilled  into  the  ces¬ 
ium  reservoir  (Fig.  3)  maintained  as  a  cold  region.  A  system  vacuum  (non- 

O 

condensibles)  of  less  than  10  torr  was  recorded  after  completing  this  opera¬ 
tion.  The  manifold  was  isolated  from  the  vacuum  system  and  Airco  reagent  grade 
argon  was  then  slowly  admitted  at  room  temperature  to  a  pressure  of  9°  torr 
(recorded  using  a  bakeable  glass  spiral  gauge)  and  the  cavity  was  pinched  off 
the  system  at  the  copper  pump-out  tubulation.  Finally  the  copper  pinch-off 
was  checked  for  leaks . 


Fig.  3*  Microwave  cavity  attached  to  the  ultra-high  vacuum  system. 

*  A  separate  oven  was  placed  over  the  cavity  during  bake-out  in  order  to  raise 
its  temperature  to  600°C  compared  with  the  450  °C  bake-out  temperature  for  the 
remainder  of  the  ultra -high  vacuum  system. 
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A  photograph  of  the  finished  assembly  including  the  coolant  tales  and 
thermocouples  is  shown  in  Fig.  k. 


Fig.  1+.  Microwave  cavity  assembly  ready  for  insertion 

into  the  3  in.  o.d.  x  4  ft  long  containment  can . 

The  circuit  for  the  measurement  of  the  cavity  resonant  frequency  is 
shown  in  Fig .  5 . 


MICROWAVE 

SOURCE 


DUAL  BEAM 

OSCILLOSCOPE  RAT|0 


RESONANT 

CAVITY 


Fig .  5 .  Microwave  circuit  for  measurement  of  cavity  resonant 
frequency . 
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Additions  to  the  circuit  described  previously ^ ^  include  a  new  wide  sweep 
K-band  microwave  generator  (H.P.696A),  a  ratiometer  (H.P.210A)  for  direct 
measurement  of  the  VSWR  and  display  on  the  oscilloscope,  new  wide  band  crystal 
detectors  (H.P.422A)  and  a  precision  attenuator  to  permit  study  of  the  shift 
of  the  cavity  resonant  frequency  (plasma  conditions)  with  variation  of  the 
incident  microwave  power  input. 

During  the  inpile  runs  on  these  cavities  with  the  new  wide  sweep  micro- 
wave  generator  we  discovered  other  resonances  as  shown  in  Fig.  6. 


Fig.  6-  Topical  signal  output  versus  frequency  sweep  (horizontal  scale  is 
22.0  to  26.5  GHz)  at  operating  conditions  (ixe-.J7xl011cm"3)  The 
upper  beam  is  the  reflected  signal  crystal  diode  output  (-V)  and 
the  lower  beam  is  the  ratiometer  output  (+V).  The  five  resonant 
peaks  are  identified  in  Table  I . 

Computations  for  both  the  TE  and  1M  modes  identified  these  resonances  as  shown 
in  Table  I.  Best  agreement  between  computed  and  measured  mode  frequencies  was 
obtained  by  assuming  d=2.287  cm  and  h=0.693  cm.  These  values  are  consistent 
since  the  originally  measured  d  =  2 .272  cm  should  increase  by  thermal  expan¬ 
sion  inpile  and  h  (0.7)  cm  could  either  decrease  or  increase  dependent  on  the 
accuracy  of  the  Kovar  cap  assembly. 
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TABLE  I .  Observed  resonances  . 


^(computed) 

f  (measured) 

Position  in 

Mode 

GHz 

(empty  cavity )GH: 

z  Fig  .6,  cm 

210 

21.44 

— 

2.0 

22.96 

22.91 

5  -2 

^020 

23.05 

23.10 

5.7 

™on 

23.85 

23.86 

6.8 

25.ll 

25.12 

9.3 

The  microwave  cavity  was  supported  by  the  end  of  the  waveguide  inside 
an  aluminum  containment  tube  as  shown  in  Ref .  1  (for  the  first  inpile  micro- 
wave  on  neon-argor).^  The  same  method  was  used  for  the  neon-argon  inpile 
run  to  be  described  here  but  in  addition  a  l/4  in.  aluminum  tube  was 
installed  so  that  it  terminated  very  close  to  the  outside  surface  of  the 
uranium  wall.  This  was  used  to  direct  a  cooling  stream  of  Ng  against  the 
uranium  wall  thereby  enabling  us  some  control  of  the  wall  temperature 
(about  80  watts  of  heat  are  liberated  in  the  1  mil  uranium  foil  at  full 
reactor  power  of  2  mW).  For  the  Ar-Cs  cavity,  a  second  Ng  coolant 
stream  was  directed  against  the  cesium  bath  in  order  to  control  separately 
the  temperature  of  the  cesium  bath.  Thermocouples  were  attached  to  the 
uranium  and  Kovar  sides  of  the  microwave  cavity  and  also  to  the  cesium  bath. 

A  photograph  of  the  operation  at  the  University  of  Michigan  research 
reactor  is  shown  in  Fig.  7.  The  microwave  generator  is  on  the  wooden  plat¬ 
form  in  the  center  of  the  reactor  bridge  and  the  waveguide  connects  to  the 
aluminum  containment  can  in  the  reactor  core.  The  ratiometer  is  at  the 
right  of  the  bridge  (above  grill)  and  the  oscilloscope  is  at  the  extreme 
right  of  the  picture.  The  aluminum  containment  can  containing  the  cavity, 
the  microwave  generator  and  the  interconnecting  plumbing  were  all  sup¬ 
ported  from  one  screw  jack  so  that  the  cavity  could  be  positioned  verti¬ 
cally  in  the  core  to  vary  the  neutron  flux  independent  of  the  reactor  power . 
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RESULTS 


A. 


Presented  in  this  section  are  the  data  from  the  second  neon-argon  micro- 
wave  cavity  and  the  data  from  the  argon-cesium  microwave  cavity. 

Electron  Density  in  Neon-Argon  Plasma 

In  the  first  run  with  this  cavity  at  low  reactor  power  the  shift  in  fre¬ 
quency  indicated  electron  densities  which  were  very  low  -  in  fact  they  were 
less  than  that  predicted  for  pure  neon.  However  as  the  run  progressed  the 
electron  density  increased  with  time  and  finally  stabilized  at  electron  den¬ 
sities  consistent  with  theory.  This  behavior  has  been  attributed  to  clean-up 
(probably  by  the  uranium)  of  a  contaminant  (possibly  0g)  which  was  liberated  at 
the  beginning  of  the  run  from  the  cavity  walls  by  the  nuclear  radiations  . 

During  subsequent  runs  there  was  an  additional  small  increase  of  electron 
density  with  time  which  may  be  due  to  a  redistribution  of  the  uranium  on  the 

foil  surface  (see  discussion  of  argon-cesium  data). 

The  data  from  the  last  run  are  shown  in  Fig .  8  where  the  electron  den¬ 
sity  is  plotted  against  the  neutron  flux.  The  solid  curves  were  computed  from 
our  reaction  kinetics  theory  and  no  adjustable  parameters  are  involved. 

Curves  were  computed  for  10"4Ar/Ne  at  90  torr  for  gas  temperatures  of  300°K 
and  600 °K.  The  only  temperature  dependence  which  we  have  built  into  the  reac¬ 
tion  rates  to  date  is  that  for  diffusion  («*  T)  and  that  for  collisional  radia¬ 
tive  recombination  («  T'9/2) .  Also  shown  for  comparison  is  a  computed  curve 

(dashed)  for  pure  neon  at  90  torr  and  300°K. 

The  experimental  data  are  shown  as  points .  With  only  nuclear  heating 
of  the  cavity,  there  was  a  gradual  increase  of  the  average  cavity^temperature 
(  =  ^Ga^  )  with  increasing  neutron  flux  from  about  320  °K  at  (j)«*10  cm  sec 
to  aboutV500°K  at  ^aHO^cm^sec  \  The  highest  electron  density  was 
1.0xl012cm“3  at  full  reactor  power  (2  mW).  The  agreement  of  experiment  with^ 
theory  is  considered  good .  The  one  experimental  point  reported  last^year^  ' 
for  the  same  gas  mixture  but  at  a  pressure  of  l40  torr  was  ng=5  -3x10  cm  at 
(|)=4xl0'1'2cm"2sec  ^  and<TQas>  =551  °K. 

At  low  electron  densities  (■*  lO^^cm  ^)  the  theory  predicts  that  the  elec¬ 
tron  density  will  decrease  with  increasing  gas  temperature  because  of  increased 
diffusion  losses  while  at  higher  electron  densities,  where  the  collisional- 
radiative  recombination  rate  is  the  predominant  loss,  the  electron  density 
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Comparison  of  experimental  and  theoretical  values  of  electron  density  versus  neutron 

flux  for  neon-argon.  The  theoretical  curves  are  for  constant  temperatures  of  300TC  and  600  K 


increases  with  the  gas  temperature.  By  increasing  the  microwave  power,  which 
heats  the  electrons  (and  thereby  the  gas),  the  electron  density  in  general 
followed  the  direction  predicted  by  theory.  However  at  high  electron  densi¬ 
ties  the  coupling  of  the  cavity  to  the  waveguide  had  become  so  small  that 
variation  of  the  power  level  of  the  incident  signal  made  only  a  small  change 

in  the  electron  density. 

The  variation  of  electron  density  with  the  cooling  gas  flow  to  the  cavity 
(and  thus  the  gas  temperature)  also  in  general  followed  the  direction  pre¬ 
dicted  by  theory  but  the  magnitude  of  the  change,  particularly  at  high  elec¬ 
tron  densities,  was  very  much  less  than  that  predicted  by  the  theory.  A  study 
is  under  way  to  include  the  temperature  dependence  of  the  other  potentially 
important  rates  in  the  reaction  kinetics  code  since  there  may  be  compensating 

processes . 

Electron  Density  in  an  Argon-Cesium  Plasma 

The  additional  parameter  that  could  be  varied  with  the  argon-cesium 
cavity  was  the  ratio  of  the  density  of  the  minor  to  major  gas  species.  This 
tube  was  filled  to  90  torr  of  argon  and  the  cesium  to  argon  atom  ratio  (Cs/Ar) 
could  be  varied  inpile  from  10”6  to  10"3  by  adjusting  the  temperature  of  the 
cesium  reservoir.  In  the  first  reactor  run  a  maximum  in  the  electron  density 
was  obtained  at  a  Cs/Ar  very  nearly  at  the  ratio  predicted  by  theory.  However, 
later  runs  showed  that  the  electron  density  was  strongly  dependent  upon  the 
temperature  of  the  cavity  walls,  and  the  wall  temperatures  had  not  been  held 
constant  during  the  first  run.  The  first  runs  were  repeated  with  care  to  main, 
tain  the  average  cavity  wall  temperature  constant .  Since  only  one  wall  (uran¬ 
ium)  of  the  cavity  was  cooled  with  gas,  it  was  not  possible  to  maintain  inde¬ 
pendent  temperature  control  of  both  walls  of  the  cavity.  The  gas  temperature 
effects  to  be  discussed  might  be  due  to  surface  phenomena  on  the  inside 
walls  of  the  cavity  and/or  collisional  processes  in  the  plasma.  Nevertheless 
we  define  the  average  gas  temperature  to  be 

y  =  1/2  ^(uranium  electrode)  +  T  (Kovar  electrode^J  (3) 

Sasav 

and  reference  to  the  "effect  of<T >  "  should  be  understood  to  include  pos- 

av 

sible  surface  reactions  . 
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Effect  of  Cs/Ar  on  Electron  Density  -  To  find  the  effect  of  Cs/Ar  on  the 
electron  density  (n  )  at  constant  gas  temperature  two  inpile  runs  were  made: 
the  first  at  a  neutron  flux  (j>=0  •72xl01^cm  2see  1  and  <^'gas)  =576  K,  and  the 
second  at  (|>=1 .44x10^  cm"  2sec  1  and  <Tgas>  =644  °K.  The  data  are  plotted  in 
Fig.  9  together  with  the  curves  predictelVfrom  the  reaction  kinetics  code  for 

the  same  conditions . 

The  electron  density  maxima  in  the  experimental  curves  occur  at  Cs/Ar 
ratios  considerably  lower  than  predicted  by  theory.  The  maxima  are  predicted 
for  Cs/Ar si.3xl0~^  while  the  experimental  maxima  are  at  an  order-of -magnitude 
smaller.  For  the  lower  experimental  curve  of  Fig.  9,  the  electron  density 
at  the  maximum  is  about  eq„al  to  that  predicted;  however  for  the  higher  experi¬ 


mental  curve  the  measured  maximum  electron  density  is  about  twice  that  pre¬ 
dicted  .  The  difference  in  electron  density  of  the  two  theoretical  curves  is 
due  mostly  to  the  difference  in  neutron  flux  and  the  large  difference  in 
electron  density  of  the  two  experimental  curves  is  evidence  of  an  anomaly  in 
the  dependence  of  electron  density  on  /T ^  • 

x  &  q  ir 


Effect  of/T  p  on  Electron  Density  -  To  establish  the  dependence  of 

the  electron  density  on  the  average  gas  temperature  (and/or  cavity  wall  tem¬ 
perature  ) ,  three  runs  were  made  at  different  constant  values  of  Cs/Ar  in 
which  the  cooling  gas  flow  rate  to  the  uranium  wall  of  the  cavity  was  varied, 
and  the  cooling  gas  flow  to  the  cesium  reservoir  adjusted  to  obtain  the  desired 
temperature  distribution  for  a  constant  Cs/Ar  ratio  in  the  cavity.  In  order 
to  span  a  significant  temperature  range  for  each  run  it  was  also  necessary  to 
make  the  runs  at  different  values  of  the  neutron  flux.  The  results  of  these 


measurements  are  shown  in  Fig  .  10  together  with  the  predicted  curves . 

In  all  three  runs  the  electron  density  increased  more  rapidly  with  tem¬ 
perature  than  that  predicted .  The  increase  in  n0  with<J?ga  \  was  greatest 
at  the  lower  cesium  concentrations  of  Cs/Ar=1.0xl0  and  5^6x10  .  Also,  at 

these  lower  cesium  concentrations  the  measured  electron  density  tended  to  be 
higher  than  the  predicted  value,  while  at  the  higher  cesium  concentration 
(Cs/Ar=10~^)  the  measured  electron  density  was  less  than  that  predicted.  Even 
at  the  lowest  cesium  concentration  (Cs/Ar=10  ^),  however,  the  measured  and 
theoretical  curves  cross  at  an  average  gas  temperature  of  4lO°K.  This  depend¬ 
ence  of  the  measured  electron  density  on  the  average  gas  temperature  (and/or 
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Fig.  9.  Comparison  of  experimental  and  theoretical  values  of  electron  density  versus  Cs/Ar  ratio  for 
two  constant  values  of  neutron  flux  and  average  gas  temperature . 


Fig .  10 .  Variation  of  electron  density  with  average  gas  temperature . 

The  theoretical  and  experimental  values  are  compared  under 
three  separate  conditions  of  constant  Cs/Ar  ratio  and 
neutron  flux. 


Fig.  11.  Experimental  values  of  electron  density  versus  reciprocal 
values  of  average  gas  temperature .  All  the  data  are  con¬ 
sistent  with  an  "activation  energy"  of  0.22  eV  represented 
by  the  slope  of  the  solid  lines . 
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cavity  wall  temperature)  is  extremely  large;  e.g.,  for  the  run  with  Cs/Ar =10  , 

an  increase  in<T  of  only  200°C(=0.017  eV)  produced  an  order-of -magnitude 

increase  in  the  SeStron  density.  In  Fig.il  the  dependence  of  the  electron 
density  on  temperature  is  shown  as  a  function  of  the  reciprocal  of  the  average 
gas  temperature,  .  All  of  the  data  from  the  runs  at  the  three  differ- 

ent  Cs/Ar  ratios  fit  t $e  exponential  function 

n  =  nQ  exp  [-£/(kT/e)j  (^) 

with  an  "activation  energy"  £=0.22  eV.  The  significance  of  £will  he  dis¬ 
cussed  later . 

The  electron  density  computed  from  the  present  model  is  already  strongly 
temperature  dependent  since  for  values  of  Cs/Ar^lO'3,  the  predominant  loss 
mechanism  for  the  cesium  ion  (CsO  is  the  3-tody  collisional-radiative  recom¬ 
bination  which  varies  as  T‘9/2.  Besides  the  T'9'2  dependence  of  this 
reaction,  the  only  other  temperature  dependence  built  into  the  present  reac¬ 
tion  kinetics  model  for  Ar-Cs  is  the  variation  of  the  diffusion  coefficients 
with  temperature .  The  temperature  dependences  of  the  other  reactions  in  the 
model  are  being  evaluated  but  individually  their  temperature  dependence  is 
expected  to  be  small  (aT1/2) .  For  this  reason  we  have  initially  sought  other 
explanations  of  this  effect.  These  are  discussed  in  the  next  section. 

DISCUSSION  OF  ANOMALOUS  TEMPERATURE 
DEPENDENCE  IN  Ar-Cs 


The  striking  influence  of  small  changes  in  the  average  cavity  tempera¬ 
ture  upon  the  electron  density  in  Ar-Cs  led  us  first  to  reexamine  our  tech¬ 
nique  of  calculating  the  cesium  density  in  equilibrium  with  the  cesium  reser¬ 
voir  -particularly  since  the  anomalous  temperature  effect  had  not  been  pre¬ 
viously  observed  in  the  non-condensable  neon-argon  system.  Additionally,  we 
considered  other  possible  processes  not  presently  included  in  the  reaction 
kinetics  equations .  What  follows  is  a  listing  with  comments  on  these  pro¬ 
cesses,  none  of  which  acting  independently  appear  likely  to  offer  an  explana¬ 
tion  of  the  observed  temperature  behavior  over  the  full  range  of  experimental 
conditions.  The  processes  discussed  are:  (a)  Dependence  of  Cs/Ar  on  <^gas^v 
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(b)  Time  to  reach  equilibrium  (c)  Delayed  surface  equilibrium  (d)  Langmuir 
surface  ionization  of  cesium  on  cavity  walls  (e)  Impurities  (f)  Inter¬ 
action  of  surface  adsorbed  cesium  atoms  with  plasma  species  (g)  Ionization 
of  cesium  by  argon  excited  states  (h)  Heteronuclear  ions  and  (i)  Influence 
of  energetic  electrons  from  Penning  ionization  process . 

(a)  Dependence  of  Cs/Ar  on<Tga4v  -  At  equilibrium  the  vapor  pressure  of  the 
cesium- in  the  cavity,  p(Cs),  is  i^ual  to  that  in  the  cesium  reservoir  and  is 
fixed  by  the  temperature  of  the  cesium  reservoir,  T(Cs)  .  However,  the  number 
density  of  cesium  atoms  in  the  cavity,  n^Cs),  is  also  a  function  of  the  gas 
temperature  in  the  cavity.  Also  the  number  density  of  argon  atoms  in  the 
cavity,  n^Ar),  is  a  function  of  <Tga^y  as  well  as  the  initial  filling  pres¬ 
sure  (p  )  since  when  the  cavity  is  heated  to  a  temperature  higher  than  the 
cesium  reservoir  temperature  some  of  the  argon  atoms  are  driven  from  the 
cavity  into  the  reservoir.  The  equation  which  takes  into  account  these 
effects  is  derived  in  Appendix  A  and  is 

Cs  eM  195  *  105  ]  (A-8> 

A?  *  P0  <V4v 


(b)  Time  to  Reach  Equilibrium  -  A  change  in  temperature  of  the  cesium  reser¬ 
voir  produces  a  change  in  vapor  pressure  in  the  reservoir.  This  leads  to  a 
flow  of  cesium  between  the  reservoir  and  cavity  to  establish  not  only  a  new 
equilibrium  concentration  in  the  gas  in  the  cavity,  but  also  a  new  equilibrium 

surface  coverage  on  the  inside  walls  of  the  cavity.  With  an  argon  gas  fill¬ 

ing  of  about  100  torr,  the  mean  free  path  of  a  cesium  atom  is  much  less  than 
the  dimensions  of  the  reservoir,  cavity  or  throat  and  therefore  diffusion 
theory  is  applicable  .  In  Appendix  B  a  rate  equation  is  derived  for  the  dif¬ 
fusion  of  cesium  from  the  reservoir  to  the  cavity  and  the  result  is  summa¬ 
rized  here.  If  n^Cs)  is  the  cesium  density  in  the  cavity  and  n2(Cs)  that 

above  the  fluid  in  the  reservoir,  we  get 


L 

[ 

I 


n^(Cs) 

n^TEsT  “ 


where  time  constant  t  is  given  by 

3(1  +  Kg)  n(Ar)  Qp  LV 

T  =  - - - - 
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(B-10) 


(B-9) 
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This  equation  assumes  that  the  gas  and  surrounding  surfaces  in  the  cavity 
approach  their  equilibrium  values  and  that  K2  atoms  of  cesium  are  deposited 
on  the  surfaces  per  unit  cesium  atom  increase  in  the  gas.  n(Ar)  is  the  argon 
number  density,  is  the  momentum  loss  cross  section  for  collisions  between 
cesium  and  argon  atoms  at  average  relative  velocity  £ •  Lis  the  length  and 
a  the  cross  sectional  area  of  the  throat  and  V  is  the  volume  of  the  cavity. 

For  typical  conditions  i  =0.675  hr  which  means  that  startlnS  from  n0 
cesium  in  the  cavity,  90%  of  the  equilibrium  concentration  would  be  reached 
in -.1-1/2  hr.  After  equilibrium  has  once  been  reached  and  conditions  changed, 
the  new  equilibrium  state  should  be  reached  within  one  hour .  Our  inpile 
experiments  support  this  conclusion. 

(c)  Delayed  Surface  Equilibrium  -  The  plot  of  the  electron  density  versus 
m  (Fig  .ll)  implied  a  volume  or  surface  reaction  with  an  activation 

energy  of  £=0.22  eV.  This  suggested  that  possibly,  contrary  to  the  assump¬ 
tions  above,  the  cesium  density  in  the  cavity  was  not  in  equilibrium  with 
that  in  the  cesium  reservoir  even  after  time  periods  of  about  1  hour,  but 
instead  that  a  temporary  additional  source  of  cesium  existed  on  the  cavity 
walls  (e.g.,  note  the  flange  on  the  uranium  wall  in  Fig.  2).  The  release 
rate  of  this  excess  cesium  would  then  depend  upon  the  cavity  wall  tempera¬ 
tures.  At  the  lower  cesium  to  argon  ratios  this  source  would  feed  addi¬ 
tional  cesium  to  the  gas  to  produce  a  cesium  to  argon  ratio  higher  than 
that  computed  from  equilibrium  with  the  cesium  reservoir  temperature.  This 
explanation  therefore  might  account  for  the  apparent  shift  of  the  maxima 
in  n  to  lower  Cs/Ar  in  Fig.  9  and  also  might  account  for  the  anomalous 
temperature  dependence .  On  the  other  hand  this  explanation  would  require  a 
means  for  the  build-up  of  excess  cesium  on  the  cavity  walls.  Such  a  build-up 
might  occur  during  the  long  down  time  (few  weeks  storage  in  the  reactor  pool) 
between  runs  inpile .  It  was  found  that  cesium  does  build  up  on  the  inside 
surface  of  the  ceramic  microwave  window  and  it  requires  about  l/2  hour  of 
inpile  operation  at  maximum  cavity  temperature  to  clear  the  window  so  that 
the  resonant  structure  can  be  seen  in  the  reflected  signal  (see  Fig.  6). 

In  order  to  check  the  above  hypothesis,  an  extended  inpile  run  (46 
hours)  was  made  at  constant  conditions  to  see  if  such  a  source  of  cesium 
could  be  exhausted .  The  conditions  chosen  were  those  for  the  maximum 


electron  density  (experimental)  for  curve  A  in  Fig .  9  so  that  depletion  of 
the  source  would  be  evidenced  by  a  drop  in  the  electron  density  as  the  Cs/Ar 
returned  from  the  high  non -equilibrium  value  to  the  lower  computed  (equilib¬ 
rium)  value  .  There  was  no  decrease  in  electron  density  for  the  48  hour  run 
and  it  was  concluded  that  no  excess  source  of  cesium  was  present.  During  ^ 
the  run  the  electron  density  actually  increased  from  1.3x10  cm  to  1.8x10 
cm"3.  This  increase  was  attributed  to  an  increase  in  gas  pressure  (fi  cion 
release  of  xenon  and  krypton)  and/or  a  redistribution  of  the  uranium-235  on 
the  walls  of  too  cavity.  According  to  the  Reactor  Operating  Staff  and  our 
neutron  flux  probe the  neutron  flux  did  not  change  during  the  run. 


(d)  Langmuir  Surface  Ionization  of  Cesium  on  Cavity  Walls  -  Because  of  the 
low  ionization  potential  of  the  cesium  atom  (Vj=3  .89  eV),  it  is  readily  ion¬ 
ized  on  a  hot,  bare  metallic  surface  particularly  when  the  work  function  is 
greater  than  3-89eV.  If  such  an  additional  source  of  ions  occurred  at  the 
walls  of  the  microwave  cavity  the  rate  of  ionization  would  be  dependent  upon 
the  temperature  of  the  walls.  The  fraction  of  arriving  cesium  atoms ^that 
are  ionized  at  a  surface  (p)  is  given  in  terms  of  Saha  Equilibrium 


1  +  2  exp  (Vj  -  $)/kT 

At  the  temperature  of  interest  in  the  microwave  studies  (400-700  K)  the  metal¬ 
lic  surface  is  partially  covered  with  cesium  atoms  adsorbed  on  the  surface . 
This  fraction  of  a  monolayer  ($)  is  itself  temperature  dependent  and,  in  fact, 
produces  most  of  the  temperature  dependence  in  Eq.(5)  at  700°K.  In  order  to 
estimate  the  possible  importance  of  surface  ionization  on  the  behavior  of  the 
microwave  cavity  it  was  decided  to  calculate  the  temperature  the  walls  must 
assume  (T')  so  that  the  surface  production  rate  of  ions  would  equal  the  volume 
production  rate  of  ions .  The  dependence  of  0  on  T  was  not  known  for  Kovar  or 
uranium  so  values  for  tungsten^  vere  used.  T'  should  be  even  higher  for 

uranium  or  Kovar.  ^ 

For  Cs/Ar=10-^  and  an  argon  pressure  of  90  torr  N(Cs)=2.9xlO  cm  and 
the  cesium  arrival  rate  R=NvA^l°l8cm"2sec"1 .  The  volume  production  rate 

*  A  surface  deposit  with  lifetime  =- 46  hours  is  possible  but  very  improbable 
since  sufficient  time  has  elapsed  to  transport  about  150  monolayers  of 
cesium  to  the  cesium  reservoir  for  a  non-equilibrium  (Cs/Ar)’  of 4x10 
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1  ^  ^ 

of  ions  equals  the  volume  ionization  rate  (S^IO  ions  cm  sec  )  times  the 
cavity  volume  (V=2.83  cm^)  and  the  surface  production  rate  equals  p  times  the 
cesium  arrival  rate  (R=10  cm"  sec  )  times  the  surface  area  (A=13 .10cm  ) . 
Equating  these  production  rates  and  solving  for  p 

SV  -3 

3  =  jfc-  =  2.2  x  10  \ 

Now  solving  Eq.(5)  for  T  gives  T=l4-50°K.  This  is  considerably  higher  than 

the  maximum  temperature  measured  so  it  was  concluded  that  Langmuir  surface 

ionization  of  cesium  was  negligible . 

The  uranium  foil  is  heated  by  nuclear  fission.  Most  of  the  fission 

energy  is  deposited  along  the  fragment  track  in  electronic  collisions  and 

it  might  be  supposed  that  very  hot  spots  would  be  produced  where  the  frag- 

( 9 ) 

ment  leaves  the  uranium  surface.  However,  investigations  have  shown  that 
most  of  the  energy  can  be  dissipated  by  the  lattice  structure  in  very  short 
times  of  order  10-11  sec.  Therefore  these  local  spots  would  not  have  suf¬ 
ficient  lifetime  and  size  to  produce  appreciable  local  Langmuir  surface 
ionization . 

(e)  Impurities  -  Unknown  impurities  are  often  invoked  as  the  cause  of 
unpredictable  plasma  behavior.  Even  so,  we  do  not  favor  such  a  cause  in  the 
present  experiments  for  the  following  reasons  . 

Stringent  conditions  of  cleanliness  were  observed  during  fabrication 
and  processing  of  the  Ar-Cs  cavity.  For  example, during  bake-out  the  cavity 
was  heated  to  600°C  and  evacuated  to  6xlO"9  torr.  The  cesium  was  carefully 
flash-distilled  into  the  cavity  reservoir  under  high  vacuum  and  then  filled 
with  Airco  reagent  grade  argon.  The  final  cavity  pinch-off  was  checked  for 
leaks  using  a  high-sensitivity  mass-spectrometer  leak  detector. 

Inpile,  fission  products  released  from  the  uranium  surface  are 
a  source  of  impurity.  The  dominant  gaseous  products  are  xenon  and  kryp¬ 
ton  which  because  of  their  nature  and  release  rate  cannot  measurably  affect 
the  electron  density  in  the  Ar-Cs  plasma  over  the  time  scale  of  the  present 
experiment .  We  can  conceive  of  other  more  deleterious  impurities  such  as 
oxygen  being  released  inpile,  possibly  arising  from  fission-fragment  uom- 
bardment  of  oxide-coated  surfaces.  However,  the  strong  gettering  action 
of  both  uranium  and  cesium  is  expected  to  minimize  the  concentration  of  such 
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impurities .  More  importantly,  our  inpile  experiment  has  yielded  reproducible 
results  from  the  outset  and  in  one  instance  the  cavity  was  operated  contin¬ 
uously  inpile  for  k6  hours  (Section  C)  with  no  significant  change  in  the  tem¬ 
perature  dependence  of  the  electron  density.  This  behavior  is  contrary  to 
that  expected  from  the  progressive  build-up  of  important  impurities  with  time. 

(f)  Interaction  of  Surface  Adsorbed  Cesium  Atoms  with  Plasma  Species  -  We 
consider  in  a  qualitative  manner  whether  a  significant  and  temperature- 
dependent  charge  generation  rate  is  likely  to  arise  via  the  interaction  of 
diffusing  plasma  species  with  the  cesium  atoms  adsorbed  on  the  walls  of  the 
Ar-Cs  cavity,  for  the  conditions  indicated  previously  in  Figs.  9  and  10.  The 
species  of  interest  are  electrons,  photons,  and  various  ionic,  metastable 
and  excited  states  of  argon  and  cesium. 

The  most  important  ion  and  electron  loss  for  this  very  collision- 
dominated  plasma  (n  zt  n.^lCr" cm  ,N  sJ 3x10  atoms  cm-5)  is  volume  recombin- 

0  *r  0 

ation.  The  atomic  Ar  ion  is  rapidly  converted  into  the  molecular  Ar^  ion 
which  promptly  decays  via  dissociative  recombination  with  an  electron.  There 
are  two  dominating  avenues  of  recombination  for  the  atomic  Cs  ion?  at  low 
Cs/Ar  values,  the  Cs+  ion  is  lost  by  colli sional-radiative  recombination  with 
2  electrons j  at  high  Cs/Ar  values,  Cs  is  converted  into  the  molecular  Cs^ 
ion  (via  a  Cs+-Cs-Ar  collision)  which  disappears  by  dissociative  recombina¬ 
tion.  We  can  therefore  consider  the  dominant  electron  and  ion  loss  paths  to 
be  via  the  products  of  recombination,  viz.,  photons,  metastable  states  and 
excited  states . 

Only  a  small  fraction  of  the  meta stable  states  will  diffuse  out  from  the 
plasma  to  the  cavity  walls.  For  example,  at  gas  pressures  *<100  torr  the  prin¬ 
cipal  loss  of  argon  metastable  states  at  the  lower  Cs/Ar  values  is  by  2  and 
3-body  collisions  with  neutral  argon  atoms  which  result  in  the  production  of 
high-energy  (<^11.5  eV)  resonance  and  non-resonance  photons.  As  the  value 
of  Cs/Ar  increases,  an  increased  fraction  of  the  argon  metastables  are  des¬ 
troyed  in  collisions  with  cesium  atoms  yielding  Cs+  ions  by  the  Penning  pro¬ 
cess,  but  as  indicated  in  the  previous  paragraph,  the  eventual  outcome  of 
such  an  ion  generation  is  likely  to  be  the  production  of  lower  excitation 
energy  («3«9  eV) .  Thus  in  either  case  the  metastable  states  are  almost 
wholly  destroyed  in  the  volume  of  the  plasma  and  few  reach  the  walls . 
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The  tentative  arguments  presented  here  suggest  to  us  that  if  an  impor¬ 
tant  charge  production  process  exists  at  the  relatively  cold  surface  of  the 
cavity,  it  probably  arises  because  of  the  interaction  of  the  surface  cesium 
atoms  with  the  copious  flux  of  photons  from  the  plasma.  Now  since  the  fis¬ 
sion  fragments  in  their  passage  through  the  gas  generate  2  argon  ions  for 
each  argon  metastable  state  (these  two  species  being  responsible  for  almost 
all  the  measured  eV/ion  pair),  the  decay  of  these  species  by  various  routes 
to  produce  photons  is  not  expected  to  yield  a  total  photon  production  rate 
very  different  from  the  volume  ion  generation  rate .  Also  if  we  assume  a  very 
large  photoelectric  yield ^  from  the  cesium-covered  surface  of  0.1  electrons/ 
photon,  then  the  electron  production  rate  at  the  surface  is  -0.1  of  the  total 
ion  generation  rate  in  the  volume  and  therefore  immediately  too  small  to  explain 
the  large  discrepancies  between  the  measured  and  computed  values  of  e  tron 
density  as  shown  in  Fig  .10.  In  addition,  the  back-diffusion  of  electrons  into 
the  surface  would  further  greatly  reduce  the  above  production  rate . 

There  is  another  reason  for  concluding  that  an  electron  source  at  the 
walls  of  the  cavity  will  not  significantly  alter  the  spatially  averaged  elec¬ 
tron  density  in  the  plasma.  The  injection  of  only  electrons  at  the  plasma 
boundaries  will  have  the  effect  of  reducing  the  ion-accelerating  (for  plasma 
positive  ions)  space-charge  potential,  thereby  reducing  the  ion  loss  rate  by 
ambipolar  diffusion.  However,  since  the  dominant  ion  loss  in  the  present  case 
is  not  by  diffusion  but  is  via  volume  recombination,  the  effect  of  changes  in 

the  ion  diffusion  rate  will  be  small. 

Finally  we  note  from  dhylor  and  Langmuir’s  work  on  the  cesium  coverage 

of  tungsten,  ^  that  over  the  temperature  range  of  interest  to  us,  viz., 
t00-700oK  where  n0  is  critically  dependent  on  the  cavity  temperature,  the 
fraction  0-n  of  a  tungsten  surface  covered  with  cesium  was  insensitive  to 
temperature  for  a  cesium  arrival  rate  101  atoms  cm  sec  (as  in  the 
present  experiment).  Thus  if  Kovar  behaves  in  a  manner  similar  to  tungsten 
and  the  adsorbed  cesium  layer  is  insensitive  to  temperature,  photoelectric 
action  at  the  surface  (even  if  considerably  underestimated  in  our  previous 
arguments)  would  be  unable  to  explain  the  observed  marked  dependence  of 
electron  density  upon  temperature . 
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(g)  Ionization  of  Cesium  by  Argon  Excited  States  -  The  equations  which  we 
use  for  computing  the  electron  density  in  the  Ar-Cs  plasma  do  not  take  into 
account  any  possible  ionization  of  cesium  by  argon  excited  states  which  are 
not  metastable  .  The  following  arguments  support  the  view  that  ionizing  col¬ 
lisions  of  argon  radiating  states  with  cesium  can  be  neglected  in  comparison 
with  the  cesium  ion  production  rate  from  argon  meta stable  states . 

Most  of  the  electrons  which  become  excited  in  argon  eventually  fall  to 

the  ground  state  via  the  lowest  (3?)^ (4s )  excited  electron  configuration 

which  consists  of  two  metastable  states  (3P2II.53  e/V,  3PQ  11*72  eV)  and  two 

radiating  states  (3P  11.62  eV,  11 .81  eV) .  This  result  follows  from  the 

considerable  imprisonment  of  resonance  radiation  and  transfer  of  excitation 

l8  "3 

which  occurs  in  argon  at  the  gas  density  2.9x10  atoms  cm  )  of  the  present 
experiment.  Thus  we  need  only  consider  here  the  influence  of  the  3P  and 
IP  radiating  states  which  present  by  far  the  greatest  concentration  of  argon 
radiating  states . 

We  now  compare  the  relative  importance  of  the  following  two  processes: 

1  IS  2 

Ar*  +  Cs - ►Cs'r  +  Ar  +  e  (0*  =  1.7  x  10"J‘P  cm  )  (l) 

Arm  +  Cs - p»Cs+  +  Ar  +  e  (h  ^  10  ^  cm2)  (2) 

m 

where  we  do  not  distinguish  between  the  two  lowest  radiating  states  (Ar*)  or 
between  the  two  metastable  states  (Arm)  because  their  respective  energy  dif¬ 
ferences  are  small.  The  ionization  cross  section  a*  for  process(l)  is  that 
computed  recently  by  Sheldon.^  ^  There  are,  however,  no  published  data 
available  for  process  (2);  the  value  of  am  shown  is  that  inferred  by 
Nolan^12^  from  his  measurements  of  ionization  growth  in  an  Ar-Cs  mixture . 

We  consider  values  of  Cs/Ar  <3xl0_1+  where  the  discrepancy  between' theory  and 
experiment  (Fig.  10)  is  greatest  and  where  we  seek  an  additional  ion  gen¬ 
eration  rate  in  order  to  account  for  the  high  measured  values  of  the  elec¬ 
tron  density.  Under  these  conditions,  the  destruction  rate  of  Arm  states 
in  the  gas  by  collisions  with  argon^1^  and  cesium  atoms  yields  a  metastable 
lifetime  t  =10"  ^  sec.  The  total  metastable  generation  rate  S  at  an  argon 

m  ,Q  -  -ir  n  1  m 

atom  density  of  2.9x10  cm"  ^  is  5x10  ?cm°sec'  ,v  half  of  which  is  directly 
due  to  fission  fragments  and  half  of  which  arises  indirectly  through  ionic 
recombination  into  metastable  levels .  Thus  the  metastable  density  is 

_  _  ,„10  -3 

n  =  St  =  5x10  cm. 

m  mm 
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We  do  not  know  the  total  generation  rate  of  the  Ar  excited  states 
but  an  upper  bound  to  S*  can  be  readily  set  by  considering  that  the  total 
energy  deposited  in  the  gas  by  the  fission  fragments  produces  exclusively 
Ar*  states  of  energy  11.7eV.  This  procedure  yields  S^-=9.5xl015cm" ^sec"1 . 
The  natural  lifetime  of  the  3^  and  IP,  states  is  expected  to  be  in  the 

-8  10  ii 

range  10  — 10  sec  and  allowing  two  orders  of  magnitude  increase  in  the 

lifetime  because  of  repeated  absorption  and  re-emission  of  the  resonance 
radiation,  we  set  t*210"°  sec.  Thus 

10 

n*  =  t*  «=;10  cm  0 


This  leads  to  a  ratio  of  Cs+  production  rate  by  excited  states  and  meta¬ 
stable  states  of 


n*  a* 


n  or 
m  m 


1-7 
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0.03. 


Since  this  ratio  is  so  small  we  conclude  that  the  neglect  of  Ar  -  Cs  col¬ 
lisions  in  our  code  calculations  is  not  the  cause  of  the  discrepancy 
between  the  theoretical  and  experimental  values  of  ng . 

(h)  Heteronuclear  Ions  -  The  ions  which  we  consider  as  important  for  con¬ 
trolling  tha  electron  density  in  the  Ar-Cs  plasma  are  Ar+,Arp,Cs+  and  Cst. 

4-  d  2 

However,  the  heteronuclear  ArCs  ion  may  also  be  present  since  such  an  ion 

(lli) 

has  been  isolated  experimentally.  Unfortunately,  the  influence  of  the 

ArCs  ion  is  not  readily  estimated  because  of  the  lack  of  data  regarding  its 
production  and  loss  rates . 

For  our  experimental  conditions,  it  appears  that  the  following  three 
processes  are  most  likely  to  lead  to  the  formation  of  ArCs+: 

Cs+  +  Ar  +  Ar  - •»  ArCs+  +  Ar  (l) 

Cs+  +  Cs  +  Ar  —  ■— ArCs+  +  Cs  (2) 

Arm  +  Cs  - »>ArCs+  +  e  (3) 

Processes  (l)  and  (2)  are  postulated  on  the  basis  of  exactly  analogous  reac- 

(l5 7  +  + 

tions  proposed  by  0skamv  "  for  the  conversion  of  He  to  HeNe  in  helium 

admixed  with  traces  of  neon.  Process  (3)  has  been  given  by  Herman  and 

(14) 

Cermak  and  is  based  on  their  measured  values  of  the  appearance  potential 
of  ArCs+. 


In  Oskam's  microwave  studies  of  helium  admixed  with  traces  of  argon, 
krypton  and  xenon,  and  neon  admixed  with  traces  of  argon  and  krypton,  at 
total  pressures  5  to  25  torr,  the  only  heteronuclear  ion  considered  to 
influence  the  plasma  electron  density  was  HeNe  ,  i.e.,  the  experimental 
results  of  the  remaining  gas  mixtures  could  be  ^adequately  accounted  for  by 
the  presence  of  the  more  usual  homonuclear  (Ne*,  Arg,  etc.)  molecular  ions^. 
However,  the  existence  of  heteronuclear  noble  gas  ions  such  is  HeNe  ,  HeKr 
and  NeAr+(in  addition  to  noble-gas  alkali  metal  ions  such  as  ArCs  )  has  been 
confirmed  by  mass  spectroscopic  studiesU6)  but  these  studies  have  been  con¬ 
ducted  at  gas  pressures  -1+xlCf3  torr.  Thus  the  stability  of  most  of  the 
heteronuclear  ions  in  the  collisional  environment  of  a  high  pressure  (100 

torr)  plasma  is  not  known.  + 

We  note  from  the  above  that  the  low-pressure  existence  of  NeAr  and  the 

absence  of  any  significant  role  by  NeAr+  in  Oskam's  studies  at  much  higher, 
pressure  suggest  that  NeAr+  is  unstable  in  a  collisional  environment .  This 
conclusion  is  supported  by  our  own  computations  of  electron  density  in  the 
Ne-Ar  system,  based  on  Ne+,  He*,  Ar+  and  Ar*  ions  only,  where  good  agreement 
is  obtained  with  the  measured  values  of  electron  density  (Fig.  8) .  Since  no 
information  is  available  on  the  stability  of  ArCs  in  the  Ar-Cs  plasma,  there 
is  the  possibility  that  the  critical  observed  dependence  of  ng  upon  gas  tem¬ 
perature  (Fig .10)  could  arise  because  ArCs+  is  stable  at  low  temperatures, 

but  dissociates  at  the  higher  temperatures. 

However,  it  should  be  noted  that  the  present  computed  values  of  ng  at 
the  lower  Cs/Ar  values  and  temperatures  of  400  to  500°K  (Fig.  10)  are  much 
less  than  the  measured  values,  and  if  ArCs+  is  introduced  into  the  existing 
computations  via  processes  (l),  (2),  or  (3),  the  result  is  a  depletion  both 
directly  and  indirectly  of  the  dominant  Cs+  ion  (process  (3)  is  presently, 
assumed  to  yield  the  Cs+  ion  only)  in  exchange  for  the  molecular  ArCs  which 
would  surely  decay  rapidly  by  dissociative  recombination  with  an  electron. 
Thus,  in  this  case,  the  ArCs+  would  act  as  a  sink  to  depress  the  calculated 
electron  density  below  its  already  low  value.  On  the  other  hand,  at  the 
high  Cs/Ar  value  of  10‘3  (Fig.  10?  a  depression  of  the  computed  curve  is 
required  to  bring  theory  into  better  agreement  with  experiment,  which  tends 
to  support  the  presence  of  an  ArCs+  ion.  Also  an  association-dissociation 
reaction  of  the  type  leading  to  a  very  temperature  dependent  molecular  ion 
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fo^tion  rate  might  reasonably  be  considered  the  cause  of  the  sensitive 

dependence  of  n  upon  temperature.  ^  .  r,  +  . 

nevertheless,  the  presence  of  a  temperature  dependent  ArCs  i 
not  alone  account  for  the  results  since  an  additional  correction  o 
existing  theory  would  be  required  at  low  values  of  Os/Ar  to  bring  e  c 
puted  values  into  closer  agreement  with  experiment. 

(1)  Influence  of  Energetic  Electrons  from  Penning  Ionization  Process  -  The 
atomic  Cs+  ions  are  generated  by  collisions  of  atomic  Os  with  argon  ( 
metastable  states: 

Arm(ll.6  eV)  +  Cs  - -Cs+(3-9  eV)  +  Ar  +  e(2  7-7  eV) . 

Thus  the  ejected  electron  can  possess  energy  €  up  to  a  maximum  of 
7.7  eY  i.e.,  the  difference  between  the  excitation  energy  of  the  arg 

metastable  state  and  the  ionization  energy  of  atomic  cesium,  w^have 

nf  10-l4cm2  for  the  above  Penning  process  and  have 
n  laree  cross  section  ol  cm 

nested  tbe  additional  ionisation  of  atomic  cesium  that  can  possibly  ar  se 
from  the  7-7  eV  electrons  during  their  time  to  reach  thermal  energies, 
now  consider  the  relaxation  time  of  these  hot  electrons. 

For  electron  densities  n^lO12™3,  the  electron  relaxation  time  r 
due  to  electron-electron  and  electron-ion  collisions  is  approximately 

0.7(T)3^2 

Tr  "  2.63  n  In  A 

6  (19) 

Where  I  is  the  electron  temperature  in  1C  and  Ais  a  factor  due  to  Spitser 
which  is  the  ratio  of  the  Debye  shielding  distance  over  an  average  criti 
impact  parameter.  quantity  A*s  in  fact  just  nine  times  the  nnmber 

electrons  in  a  Debye  sphere  1  (  and  is  given  by 

v  m  3/2 

*  (12  np'  ^ 

A=  nel/2 

For  values  of  T=800°K  and  ne=1012cm-3,  mA=3-3-  ®>is  yields  a  value  for 
T  of  3.  8x10-9  sec  and  thus  we  note  that  because  of  the  strong  influence  o 

coulomb  forces,  the  hot  electrons  are  rapidly  thermalised  ^^"se 
electrons.  The  effect  of  hot  electron  -  atom  collisions  will  be  to 
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Tr  still  further  by  a  snail  amount  -  the  decrease  being  small  because  only 
elastic  collisions  are  permissible  with  argon  for  €-<11.5  eV. 

We  estimate,  very  approximately,  the  number  f^  of  Cs  atoms  ionized  by 
a  hot  electron  by  using 


where  n„  is  the  density  of  Cs  atoms,  a.  is  some  average  ionization  cross 
Cs 

section  for  e-Cs  collisions  where  the  electron  energy  degrades  from  -o  7  eV 
to  below  the  ionization  potential  of  cesium  in  times  —  Tr,  and  vg  is  taken 
to  be  the  velocity  of  a  5  eV  electron.  For  a  value  of  Cs/Ar  =10  ,  nCg= 
2.9xlOll*'cm”^  and  we  assume  a  value  for  a±  of  10*15cm2.  Thus  with  ve= 

1.3x10® cm  sec"1  and  Tr=l  .8xl0"^sec,  we  obtain  f^O.07. 

On  the  basis  of  this  result,  we  have  therefore  tentatively  concluded 
that  there  exists  no  unusually  high  tail  to  the  Maxwellian  distribution  of 
the  plasma  electrons  and  that  no  large  error  is  made  by  neglecting  the  ener¬ 


getic  electrons  from  the  Penning  process,  particularly  for  values  of 
Cs/Ar  :*10’\ 
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APPENDIX  A 


COMPUTATION  OF  CS/AR  FOR  MICROWAVE  CAVITY 


For  the  argon  pressures  of  interest  in  these  studies  the  mean  fi,  ee 
path  of  an  atom  is  much  less  than  the  dimensions  of  the  cavity  so  we  use 
the  ideal  gas  law  P=nkT.  At  the  conditions  of  seal-off  the  gas  pressure  is 
p0  and  the  temperature  is  TQ  (-^  300°K).  Therefore  we  have 

Po(Ar)  =  nQ(Ar)kTo  , 

At  operating  conditions,  the  cavity  of  volume  ^  is  at  a  temperature  ^ 

N  )  and  the  cesium  reservoir  of  volume  V  is  at  a  temperature  Tg. 

'  NSasav 


p^Vi 


Resonant  Cavity 


P2'T2,V2 


Cesium  Reservoir 


From  the  conservation  of  the  total  number  of  argon  atoms  and  assuming  zero 
thermal  expansion  of  the  volumes 

p 

n1(Ar)V1  +  n2(Ar)V2  =  n0(Ar)(V1  +  Vg)  =  (VV  (A-2 


Now  p^(Ar)=Pg(Ar)  so  from  the  ideal  gas  law 

n2(Ar)  =  \  ni(Ar) 


(A-3) 


Tr 


Substituting  A-3  into  A-2  gives  the  number  density  of  argon  in  volume  1 


nl(Ar)  = 


V  + 

1  2 

1  — tt“” 

0  vi  +  T  '2 
i  ■l2 


(A-4) 


At  steady  state  the  vapor  pressure  of  cesium  in  is  the  same  as  that  in  Vg 
and  is  fixed  by  the  temperature  (Tg)  of  the  liquid  cesium  in  Tg. 

5l(0»)  =P2(Cs)  =  f(Ta).  (A-5) 
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Also  from  the  ideal  gas  law 


p^(Cs)  =  n1(Cs)kT1 


(A-6) 


Solving  A-6  for  n^Cs)  and  dividing  by  n^A r)  from  (A-4)  gives  the  cesium-to- 
argon  ratio  in 


Cs 

Ax 


n^Cs)  i  P1(Cs) 


■I  P 

I  it,  IK  + 


T„ 


V  + 

1  2 


(A-7) 


For  the  argon-cesium  cavity  V1=3.39  cm3,  Vg=1.8l  cm3  and  the  volume  fractions 
are  therefore  V^/(V^+Vg)=0 .65  and  V^/(V^+V^)=0 .3^ .  For  T_=300°K  Eq.(A-7) 
becomes 


Cs 

p(Cs) 

195  x  105 

Ar  = 

*0 

.^ga^av  TCs. 

(A-6) 


t 
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APPENDIX  B 


DIFFUSION  OF  CESIUM  IN  AEGON 


In  the  operation  of  the  argon-cesium  microwave  cavity,  tne  time  to 
reach  equilibrium  (i.e.  steady-state)  conditions  depends  primarily  upon  the 
diffusion  rate  of  cesium  (in  argon)  between  the  cavity  and  the  cesium  reser¬ 
voir  .  A  change  of  vapor  pressure  of  cesium  in  the  reservoir  leads  to  a 
flow  of  cesium  between  the  reservoir  and  cavity  to  establish  not  only  a  new 
equilibrium  concentration  in  the  gas  in  the  cavity  but  also  a  new  equilibrium 
surface  coverage  on  the  walls  of  the  cavity.  In  the  following  a  diffusion 
rate  is  computed  for  an  assumed  non-equilibrium  condition  and  then  a  time  to 
reach  equilibrium  is  computed  for  a  case  which  starts  considerably  removed 
from  equilibrium. 

The  equation  for  diffusion  along  a  tube  has  been  derived  previously' 

and  is  - 

T  _  1  JLL.  (B-l) 

Cs  3  Vd  dz 

where  I  is  the  diffusion  current  passing  through  the  cross  sectional  area  a 
of  the  tube,  u  is  the  average  relative  velocity  for  collisions  of  argon  and 
cesium  atoms,  is  the  momentum  loss  cross  section,  is  the  concen¬ 

tration  of  argon  atoms  and  dn(Cs)/dz  is  the  gradient  of  the  concentration  of 

cesium  atoms  (n^)  along  the  axis  of  the  tube  (z). 

u  and  n^)  are  actually  functions  of  z  but  for  this  computation  we 

assume  them  independent  of  z  and  also  approximate  the  gradient  as 

dnca  _  [na(Cs)  -  n!(Cs)  1  (B-2) 

where  n2(Cs)  is  the  cesium  density  in  the  reservoir,  n^Cs)  is  the  cesium 
density  in  the  cavity  and  L  is  the  length  of  the  tube  joining  the  reservoir 
to  the  cavity.  The  pertinent  dimensions  of  the  device  are  given  in  Fig.  B-l 
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Total  interior  cavity  area: 

S  -  2  [f  x  (2.272)2]  +  [it  x  2.272  x  O.7J  =  8.10  +  5  .0  =  13.10  cm2 


Fig.  B-l .  Dimensions  of  Ar-Cs  microwave  cavity. 

We  select  as  topical  the  following  temperatures:  T(ur)  =  ^20°C, 
T(k)=300°C  and  T(Cs)=l46°C  and  therefore  <Tga^v.=(U20+300)/2=360oC=633°K 
and  the  average  gas  temperature  in  the  connecting  tube  is<T>=(  36CH-lh-6)/2= 
253°C=526°K. 


The  average  relative  velocity  is 


(1) 


u  = 


'8^%  +  mCs) 


1 1/2 


mAr  mCs 


(B-3) 
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and  substituting  =526 °K  and  the  masses  of  argon  (m^)  an^  cesium  (m^,  )  we 
obtain  u=5.98xl0^cm  sec  ^ .  The  maximum  diffusion  current  occurs  when  n^(Cs)=0 
that  is,  immediately  after  the  reservoir  is  brought  to  temperature  and  before 
any  cesium  has  reached  the  cavity.  At  To(Cs)=l460C=4l9°K,  p_(Cs)=7 *0x10  1  torr 

t]  I  ^ 

and  n2(Cs)=p_(Cs^^.Tg(Cs)=l .61x10  cm"  .  Substituting  these  values  and 

Q  =5  .72x10" 1  cm^  into  Eqs.(B-3)  and  (B-l)  we  obtain  for  the  maximum  diffu- 

^  12  —1  15  —1 

sion  current  ICg (max) =2 .48x10  atoms  sec”  =8.95x10  }  atoms  hr 

When  the  partial  pressure  of  cesium  in  the  cavity  equals  that  in  the 

reservoir  the  total  number  of  cerium  atoms  in  the  cavity  is  N.  (Cs)=n..  (Cs)V,  = 

-i  Ji  -)  t.  XXX 

n2(Cs)V1T2/T1=l .61x10  x3 .39x419/526=4 -34x1 r-  atoms  of  cesium.  At  the  maxi¬ 
mum  rate  the  time  to  saturate  the  gas  would  only  be  t=N^(Cs)/l^s(max)=0 .0484  hr 

but,  of  course,  the  driving  force  tends  to  zero  as  p^(Cs) - ^P2(Cs).  The 

diffusion  rate  decreases  as  equilibrium  is  approached  but  the  adsorbed  layer 
of  cesium  on  the  inside  surfaces  of  the  cavity  walls  must  also  come  to  equi¬ 
librium  with  the  partial  pressure  of  cesium  in  the  gas . 

8  2 

For  a  cesium  atom  diameter  of  5 *3^x10"  cm,  the  projected  area  (jt/4  D  )  is 
-15  2 

2.2x10  cm  per  atom.  From  the  inverse  of  this  area,  the  atoms  in  one  mono- 


layer  of  adsorbed  cesium  is  4  ^xlO^atoms  cm"^, 


The  cesium  arrival  rate  at 


the  wall  surfaces  srfn^(Cs)  x  v(Cs)  lO^cm’^xlO^cm  sec  ^=10^  atoms  cm  ^sec  1. 
At  this  arrival  rate  and  temperature  (526°K)  the  fraction  of  the  surface 

(a) 

covered  (assuming  surface  behaves  like  tungsten)  is  ^>0.9.  Hie  inside  sur- 

2 

face  area  of  the  cavity  is  13 .10  cm  so  the  total  number  of  atoms  adsorbed 

for  0.9  of  a  monolayer  is  Nc(Cs)=0.9x4 .5x10^x13 *1=5  .3x10^  atoms  which  is 

o  14 

about  10  times  the  atoms  contained  in  the  gas  (N^(Cs)=4 .34x10  atoms). 

Overall  Rate  Equation 

An  approximate  overall  rate  equation  can  be  derived  as  follows .  From 
Eqs .  (B-l)  and  (B-2)  the  diffusion  current  through  the  throat  can  be  written 


ICs  =  4  [n2(Cs)  -  n1(Cs) j 


(3-4) 


where 


1  u  a 
3  n(Ar)Sl  L 


(B-5) 
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The  cesium  diffusing  through  the  throat  builds  up  the  concentration  in  the 
gas  in  the  cavity  and  also  deposits  on  the  walls  j  thus 


dn^(Cs) 


dns(Cs) 


(B-6) 


We  assume  the  cesium  transport  rate  across  the  cavity  to  the  walls  is  large 
compared  to  the  transport  rate  through  the  throat  and  the  cesium  deposition 
rate  on  the  surface  is  proportional  to  the  cesium  concentration  in  the  gas 
n-^Cs).  Therefore  we  set 


dng(Cs) 


dn^(Cs) 


(B-7 ) 


and  combining  Eqs  .  (B-4),  (B-6)  and  B-7) 


where 


dn  (Cs) 

—at-  *  t  h(Cs)  -  ni(Cs)] 


3(1  ♦  k2)  n(Ar)  «d  L  V 


(b-8) 


(B-9) 


Integrating  Eq.(B-8)  with  ng(Cs)  independent  of  time 

ni(Cs)  .  -  t/T 

n_tCs)  "  e 


(B-10) 


In  this  model  the  gas  and  surface  coverage  approach  their  equilibrium  values 
together  with  kg  atoms  to  the  surface  for  each  atom  to  the  gas .  t  is  the 
e-folding  time  as  equilibrium  is  approached . 

Setting  kg  =  10  cesium  atoms  on  the  surface  per  cesium  atom  in  the  gas 
and  substituting  into  Eq.  (B-9)  we  obtain 

t  *  3  x  11  x  (2.9  x  1018)  x(5 .72  x  10"lk)  x  (0.95)  x  (3-39) 


(5.98  x  10  7  x  (0.122) 
t  =  2420  sec  =  O.675  hr. 


Starting  with  no  cesium  in  the  reservoir  90  J0  of  the  equilibrium  concentra¬ 
tion  should  be  reached  after  1-1/2  hr.  After  equilibrium  has  once  been 
obtained  and  conditions  are  changed  for  a  new  run,  the  new  state  of  equi¬ 
librium  should  be  reached  within  one  hour. 
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SECTION  B 


ELECTRON  TRANSPORT  IN  A  DIFFUSION  CONTROLLED  PLASMA: 

DESCRIPTION  OF  COMPUTER  PROGRAM  AND  PARAMETRIC  STUDIES 

ABSTRACT 

The  computer  program  is  described  for  the  diffusion  model  of  the  electron 
transport  problem  in  which  the  plasma  is  generated  by  fission  fragments  and 
decays  by  ambipolar  diffusion  to  the  electrodes .  In  the  temperature  range  of 
thermionic  energy  conversion  the  upper  limit  of  applicability  of  this  model, 
where  volume  charge  loss  becomes  appreciable,  is  shown  to  be  at  an  electron  den¬ 
sity  of  2xl0'*'^cm"^  .  Parametric  studies  on  the  Richardson  currenc,  ion  genera¬ 
tion  rate  and  electrode  spacing  showed  that  for  diode  currents  of  interest  in  a 
practical  energy  converter  (l«*J-10  amps  cm'2)  the  necessary  electron  density 
was  beyond  the  limit  of  applicability  of  this  theory. 

I .  OBJECT 

To  program  the  equations  of  the  analytic  solution  of  the  diffusion  model 
of  the  electron  transport  problem  for  computation  on  the  IIM  7°9^  computer  and 
to  perform  sufficient  computations  to  ascertain  the  applicability  of  this  model 
to  our  noble  gas  thermionic  energy  converter  studies. 

II .  CONCLUSIONS 

1.  The  basic  hypothesis  for  this  model  of  charge  loss  by  ambipolar  diffusion 
is  applicable  for  low  electron  densities  and  small  spacings  (e.g.  ng  c 
2xl0^2cm"^  for  <1=0.1  cm)  where  volume  recombination  losses  are  negligible. 

2.  For  diode  currents  of  interest  in  practical  thermionic  energy  converters 
(l-sJ**20  amps  cm'2),  the  predicted  electron  densities  exceed  the  limit  of 
applicability  of  this  model  and  enter  the  range  of  appreciable  volume  charge 

loss . 

3.  A  solution  of  the  electron  transport  problem  is  needed  which  also  includes 
volume  loss  of  the  charges  wherein  a  higher  electron  density  is  available 
for  electron  transport . 
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4.  This  diffusion  model  with  a  uniform  volume  ionization  rate  from  fission 
fragments  predicts  considerable  variation  in  the  electron  density  between 
the  emitter  and  collector;  typically  ne(max)/ne(min)  -  40. 

5 .  The  calculated  diode  current  was  very  sensitive  to  the  ion  generation  rate 
and  diode  spacing  is  expected.  For  the  plasma  limited  conditions  studied, 
the  diode  current  was  relatively  insensitive  to  the  emitter  Richardson 
current . 

Ill .  INTRODUCTION 

The  output  power  from  a  thermionic  energy  converter  is  strongly  dependent 
on  the  electron  transport  properties  of  the  inter electrode  plasma .  The  con¬ 
ductivity  of  the  plasma  is  proportional  to  the  electron  number  density  so 
that  in  our  recent  studies  on  noble  gas  plasmas  generated  by  fission  frag¬ 
ments^"^  we  have  concentrated  on  maximizing  the  electron  number  density  ^b^- 
maximizing  the  ion  generation  rate^  and  minimizing  the  ion  loss  rate. 

These  studies  were  made  in  the  absence  of  thermionic  electrons .  With  high 
current  densities  of  thermionic  electrons  in  a  practical  energy  converter 
the  electron  density  distribution  across  the  plasma  is  itself  influenced  by 
the  net  current  across  the  plasma .  This  implicit  dependence  of  current  on 
electron  density  distribution  makes  the  voltage  characteristic  for  a  therm¬ 
ionic  converter  very  difficult  to  compute . 

In  an  earlier  attempt to  compute  the  I-V  characteristic  for  our  noble 
gas  thermionic  converter  at  low  current  densities,  the  transport  properties 
of  the  plasma  were  ignored  and  the  I-V  characteristic  was  computed  from  con¬ 
siderations  of  just  the  random  ion  and  electron  currents  at  the  plasma  edge- 
sheath  regions .  This  theory  assumed  a  uniform  elec  -on  density  di  stribution 
across  the  plasma  and  the  computed  I-V  curves  agreed  with  the  shape  of  the 
experimental  curves  at  low  currents  but  not  in  the  region  approaching  satura¬ 
tion  of  the  current  density . 

Recently  a  diffusion  theory  was  formulated  by  Rees  ^  for  the  transport 
of  thermionic  electrons  through  a  plasma  contained  by  parallel  plane  elec¬ 
trodes.  The  current -voltage  characteristic  was  derived  for  a  plasma  genei- 
ated  by  fission  fragments  and  which  decays  by  ambipolar  diffusion  to  the 
walls  of  the  diode .  Volume  charge  loss  was  not  included .  This  report  describes 
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the  computer  program  which  solves  the  analytic  expressions  for  the  current- 
voltage  characteristic .  Some  parametric  studies  using  this  code  are  dis¬ 
cussed  . 

IMs  diffusion  model  is  applicable  to  a  plasma  diode  operating  at  low 
values  of  pressure  and  gap  separation*where  the  electron  density  is  moderate 
(n  «  10  cm-'3).  For  the  conditions  of  a  practical  thermionic  energy  converter 

*  1  q 

(n 10  3)  the  volume  charge  loss  is  expected  to  be  appreciable  or  even  to 
dominate  in  which  case  the  diffusion  model  would  not  be  applicable. 


IV.  ANALYTIC  SOLUTION  TO  TRANSPORT  EQUATIONS  FOR 
A  DIFFUSION  CONTROLLED  PLASMA 


Transport  Equations  for  a  Collision-Dominated  Plasma 

The  motion  of  positive  ions  (of  various  species,  i)  and  electrons  at  a 
point  in  a  gas  may  be  expressed  in  terms  of  the  concentration  gradients  and 
electric  field  at  that  point  (in  which  they  are  at  equilibrium)  as: 


P "  =  n-  v"  =  -  D“Vn“  -  |i"  E  n 


ni  vi  = 


T7  +  +  _  + 

-  \  Vni  +  ^  E  \ 


^net  =  r  i 


<(l) 

(2) 

(3) 


where  I  represents  particle  current  density;  n  is  particle  density;  v  is  the 
average  drift  velocity;  D,  the  diffusion  coefficient;  n,  the  mobility;  E,  the 
electric  field;  and  Ciet  is  the  total  resultant  current.  Considerable  simpli¬ 
fication  is  obtained  with  the  assumptions  of:  X 

(1)  a  predominant  (or  an  effective)  ion  species,  n+= 

(2)  a  quasi-neutral  plasma,  (n+-n")«:  n"  so  n+=*n"=n, 

(Consequently,  from  Poisson's  equation  dE/dx=(e/£>o)(n""-  2^n^)ai0); 

(3)  diffusion  and  mobility  coefficients  which  are  independent  of  the  tem¬ 
perature  variations  across  the  interelectrode  gap  but  which  do  vary 
with  the  average  temperature  of  the  respective  particles;  and 


*The  limitation  of  this  theory  to  small  gap  spacing  means  that  these 
results  do  not  apply  to  the  microwave  cavities  used  in  the  experiments 
of  Section  A. 


(4)  parallel  plane  geometry  (V - >d/dx)  where  <  Ag«=*=d  and 

is  the  Debye  length  for  electrons ,  is  the  electroa-neutrnl  atom 
mean  free  path  and  d  is  the  interelectrode  spacing , 

With  these  assumptions  the  particle  transport  equations  simplify  to 


dn 

dx 


00 


where  D  =(D+n“+D“n+)/(|i++ti“)  is  the  ambipolar  diffusion  coefficient, 
s. 

When  volume  charge  loss  is  included,  the  continuity  equations  for  the 

charges  can  be  expressed  in  terms  of  the  ion  generation  rate  (s(r  )  )  and  an 

effective  volume  recombination  coefficient  CL  (N  n)  as 

°> 


^d7"  =  lfx  =  SW  "  a  (nM)2  ^5) 

Hie  analytic  solution  to  Eqs .  (4)  and  (5)  is  not  known  when  the  last 
term  is  not  negligible .  Although  the  model  will  be  somewhat  restrictive  it 
was  desired  to  obtain  the  analytic  solution  for  the  simpler  case  of  diffusion 
loss  only.  Consequently  with  the  further  assumptions  of 

ip  2 

(5)  negligible  volume  charge  loss  (a  n2««  (D~/ A  )n)  and 

(6)  uniform  ion  generation  rate,  S=S=  /  S(”r  )d  r/V, 

-'V 

equation  (5)  reduces  to 


4^=  4-  = s-  <«> 

dx  dx 

Substituting  T  from  (6)  into  (4)  and  integrating  gives  the  expression  for 

(3) 

the  plasma  density  across  the  plasma' 


n(x)  =  n(x 
where  -  Sx. 


(7) 


Boundary  Conditions  -  Sheath  Equations 

In  an  operating  thermionic  converter,  the  plasma  is  bounded  by  sheaths 
at  the  electrode  surfaces  which  are  negligible  in  thickness  (unignited  mode) 
compared  to  the  distance  between  the  electrodes .  The  potential  diagram  in 
Pig.  1  shows  the  thickness  of  the  sheaths  greatly  exaggerated.  Hie  electron 
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Fig.  1.  Schematic  of  potential  distribution  across  diode. 

and  ion  currents  to  and  from  the  plasma  across  the  sneaths  provide  the  bound¬ 
ary  conditions  to  complete  the  solution  of  Eq.(T)  for  n(x)  and  Eq.(l)  for  E. 
For  a  given  n<*,current  (rnet=Jnet(A/cm2)/e)  the  voltage  drop  across  the 
plasma  (Vp=- J  E  dx)  can  be  computed  and  combined  with  the  sheath  potentials 
and  electrode°work  functions  to  give  the  voltage  drop  across  the  diode .  The 
assumptions  made  here  are 

(7)  collisionless  sheaths; 

(8)  Maxwellian  distribution  of  ion  and  electron  velocities  in  the  plasma 
where 

(9)  the  ion  temperature  (T+)  an!  electron  temperature  (T‘)  of  the  plasma  at 
the  edges  near  the  sheaths  are  set  equal  to  the  corresponding  electrode 
temperature  (Tg  or  T^)  and 

(10)  the  ion,  electron,  and  gas  temperatures  in  the  main  bulk  of  the  plasma 
are  set  equal  to  the  mean  electrode  temperature  (Tp=(TE+Tc)/2); 

(11)  the  thermionic  Richardson  current  from  the  emitter  is  always  sufficiently 
greater  than  the  random  current  from  the  plasma  at  the  emitter  that  the 
emitter  sheath  is  electron-rich  and  the  emitter  sheath  potential  always 
retarding  to  ion  fow  from  the  plasma;  and 
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(12)  the  thermionic  Richardson  current  from  the  collector  Is  negligible  and 
the  collector  sheath  potential  will  depend  upon  the  magnitude  of  the  ne 
current.  (This  requires  formulations  of  the  collector  sheath  equations 
both  for  ion-retarding  to  the  plasma  and  also  electron-retarding  to  the 

plasma .) 

Analytic  Solutions  .  „ 

- -  n  x,  Vn+Vt  +hp  -ion-rich  and  electron-rich  collector 

The  analytic  solutions  for  both  the  ion-ricn  unu 

sheaths  have  been  presented  elsewhere  .(3)  These  equations,  however,  were  no 

in  a  form  that  could  be  programmed  directly  for  solution  on  the  computer. 

equations  were  reformulated  and  are  presented  in  the  Problem  Statement  in 

Appendix  A. 

The  measurable  potential  difference  between  the  electrodes  of  a  therm¬ 
ionic  converter  is  given  in  terns  of  the  difference  in  work  functions  of  the 
emitter  (»„)  and  collector  ♦„)  and  the  internal  diode  voltage  (V)  as 

Diode  Voltage  =  (*c  -  AE)  +  v 

The  electrode  work  functions  are  independent  of  the  plasma  and  sheath  poten¬ 
tials  so  we  set  V*E  aad  aalculata  1  VerSUS  V‘  **  ieSUllm  ^  charaCt"‘ 
istic  can  readily  be  corrected  for  a  work  function  difference  by  shifting  the 

curve  along  the  V-axis  by  • 

Recombination  limit  and  AppH --hi  1  it;,  of  Diffusion  Model 
When  the  volume  recombination  loss  of  ions  is  included  in  Eq.(5)  the 
analytic  solution  to  Eqs.«  and  (5)  cannot  be  found.  Nevertheless  the  error 
introduced  by  neglecting  this  term  can  be  estimated  from  the  diffusion  solu¬ 
tion  for  the  case  where  the  error  is  small  and  the  solution  would  be  only 
slightly  perturbed.  If  we  set  ah  upper  limit  on  the  error,  an  upper  hound 
can  be  computed  for  the  electron  number  density  for  which  the  diffusion  mode 

is  applicable. 

We  arbitrarily  set  the  limit  of  validity  of  the  diffusion  model  in  terms 
of  the  ion  diffusion  loss  (Lj)  and  ion  recombination  loss  (I*)  as  Ig/I*  <  £ 
where  C  is  a  small  fraction  like  0.1.  In  the  diffusion  model  L„  is  equal  to 

the  ion  generation  rate  (Sd)  so  we  have 


h .  j:°  n<; 

Sd  Sd 


(x)2dx 


<  £ 


(8) 
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Substituting  n(x)  from  Eqs.(7)  and  (8)  gives 


Also, using  Eq.(7),  the  average  electron  (or  ion)  density  across  the  gap  is 


and  L^/Sd  can  be  expressed  in  terms  of  n  from  Eqs  .  (8)  and  (12)  as 


-sr  -  s  I”2  +  Kl<  e 

(13) 

where 

0  V*  C  Ci'  C  2d2 

K  .  -L_  +  1  2  +  2  . 

(lit- 

30  3  12 

Finally  the  range  on  the  average  electron  density  (n.)  where  the  volume  recom¬ 
bination  loss  can  be  neglected  in  comparison  to  the  diffusion  loss  is 


Hie  sample  problem  in  reference  (5)  was  chosen  for  a  noble  gas  therm¬ 
ionic  converter  operating  in  the  reduced  neutron  flux  of  a  research  reactor 
30  .2  -1 

■(♦=1.0x10  Jcm  sec  )  at  selected  conditions  which  are  not  necessarily  opti¬ 
mum.  For  this  problem  S=1  .lxlO^cnT^sec"'*'  and  K=-4 ,4lxl0^cm"^ .  From  the 

f  -3  N  o  p 

reaction  kinetics  studies  for  neon-argon' 3  at  these  conditions  a'=(S-D  n/A  )/n  = 
_q  _i  a 

2.8x10  cm  sec  which  is  near  the  minimum  value.  Substituting  these  values  in 

(15)  with  £=0.1  we  find  n<  2xl0^cm"”^  and  conclude  that  this  diffusion  model 
may  be  useful  in  the  analysis  of  the  forthcoming  inpile  transport  experiment 
but  another  model  which  includes '  volume  recombination  loss  of  the  charges  will 
be  needed  to  predict  the  current -voltage  characteristics  for  a  practical  con¬ 
verter  (n/SlO^cm"^) . 
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V .  COMPUTER  PROGRAM 


This  program  computes  the  current -voltage  characteristic  from  the 
analytic  solution  to  the  diffusion  model  of  the  transport  equations.  The 
equations  for  the  analytic  solution  are  presented  in  the  problem  statement 
in  Appendix  A  and  the  flow  of  the  program  is  essentially  the  reverse  of  the 
order  of  presentation  of  the  equatioi  -  Selected  values  of  the  net  diode 
current  are  input  and  the  diode  voltage  (for  *E=«C)  is  output.  Other  input 
variables  are  listed  in  Table  A-l. 

A  listing  of  the  program  is  given  in  Table  B-l  of  Appendix  B.  The  code 
output  from  the  trial  problem  is  given  in  Table  B-2  and  a  sample  input  data 
sheet  is  shown  in  Table  B-3* 

VI .  PARAMETER  STUDIES 


After  checkout  of  the  code,  parameter  studies  were  initiated  to  investi¬ 
gate  the  influence  on  the  I-V  characteristic  of  the  magnitudes  of  the  elec¬ 
trode  spacing  (d),  the  ion  generation  rate  (s)  and  the  Richardson  current 
density  (jR)  .  The  results  of  these  runs  are  summarized  in  this  section .  A 
summary  of  the  results  is  presented  in  Table  I . 

I-V  Characteristic  for  Trial  Problem 

The  I-V  characteristic  for  the  trial  problem  (see  reference  (3)  and 
Table  A-l)  for  d-0.1  cm,  S-l.lxlO^cm^sec"1  and  JR=2.0A/cm  is  shown  in 
Fig.  2.  The  absolute  value  of  the  net  current  density  is  plotted  on  the  log 
scale  on  the  ordinate  .  The  code  will  also  calculate  the  voltage  for  nega¬ 
tive  values  of  the  current.  The  curve  is  in  agreement  with  the  hand  calcu¬ 
lated  values  of  reference  (3).  The  electron  density  distribution  across 
the  gap  is  shown  in  the  insert  of  Fig.  2  with  nQ=2 .10x10  cm  and  n^  ^ 
0.328xl01:Lcm'3.  The  average  electron  density  from  Eq.(l2)  is  0.87x10  cm 
which  is  considerably  higher  than  the  values  at  the  plasma-sheath  boundaries. 
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10.1-11 

12 .1- 14 .1 

15.1- 1?  .1 
18.1 

25.1 

24.1 
23, .1 

22.1 
28.1 
29-1 

C)  Effect  of  d 


0.1 

100.0 

0.01 

10.0 

0.05 

10.0 

0.10 

10.0 

0.50 

10.0 

1.0 

10.0 

0 .262x10 

0.385 

0.480 

0.352x10 

0.525 

0  .220x10 

O.376 
1.12 
1.69 


0 .262x10" 
0.352 
0.414 
0.462 

1 .l4xl0lc 

1.36 

1.54 

1.69 

3.07 

3.76 


0.576x10 
0.403  1.25 
0 .766  1 .69 
3.76  3.07 

7.53  3.73 


0 .115x10 


0.225x10 


0.991x10- 

0.995 

1.01 

1.03 


0 .115x10' 
0.225 
0.335 
0.445 

0 .421x10- 
0.624 
0.825 
1.03 
5.00 
9.95 

0 .194x10 
0.539 
1.03 
5.04 

1 .OlxlO1 


1 .81x10 

1.89 

7 .52x10' 

7 .65 
8.02 
8.29 


0.94x10 

1.81 

2.64 

3-48 

3-52xlO] 
5.16 
6 .69 
8.29 

3-91x10; 

7.58x10' 


0.455x10 

2.65 

8.29 

1 .93X101 
7 .04x10 


Effect  of  Richardson  Current  (JR) 

The  effect  of  varying  just  the  Richardson  thermionic  electron  current 
from  the  emitter  (JR)  on  the  saturation  current  density  (Jgat)  is  shown^in 
Table  I-A.  for  a  diode  spacing  of  d=0.1  cm  and  three  values  of  S=l.lxlO  , 
2.2xl016  and  1.0xl017cm"3sec_1.  The  effect  is  very  small  since  Jgat  is  fixed 
primarily  by  the  random  electron  current  density  at  the  collector  (providing 
random  electron  current  density  at  the  collector).  From  Eq.(A-36)  JBa^- 
q(v-)n^ax=(°.7l4.6xlO“12)nd.  The  random  electron  current  depends  strongly  on 
the  ion  generation  rate  (s)  and  this  is  reflected  in  the  increase  of  Jsa-t  as 
S  is  increased  from  l.lxlO16  to  l.OxlO^cm^sec"1 .  The  electron  densities  at 


the  plasma 

value  <n> 

-3  -1 ' 

cm  sec 


sheath  interfaces  are  recorded  in  Table  I  as  well  as  the  average^ 
calculated  from  Eq.(l2).  <n>  varies  from-^1012cm  3  at  S=l.lxlO 
to -■<■*10'*" 3  for  S=1  .OxlO'*"7cm"3sec"'*'  so  that  from  the  criterion  of  appli¬ 


cability  in  (15),  volume  recombination  loss  of  ions  would  be  appreciable  above 
S^^xlO  cm  ''sec 


Effect  of  Ion  Generation  Rate  (s) 

The  effect  of  the  ion  generation  rate  (S)  on  J^at  and  n^is  shown^further 
in  part  B  of  Table  I  where  S  is  varied  from  1 .1x10  to  1x10  cm  sec  . 

The  effect  of  the  ion  generation  rate  (s)  on  the  shape  of  the  I-V  char¬ 
acteristic  is  shown  in  Pig .  3  •  To  obtain  J>1  A  cm-2  for  d=0.1  cm  would 
require  S'°10'^cm"3sec~"*'  but  as  shown  in  Table  I  the  electron  density  for 
these  conditions  greatly  exceeds  the  limit  of  applicability  of  the  diffusion 

model,  Eq.(l5)  • 

The  open  circuit  voltage  (V  at  J=0)  increases  (jVl  decreases)  as  the  ion 
generation  rate  increases.  This  is  a  consequence  of  the  increased  random  elec¬ 
tron  current  from  the  plasma  to  the  emitter  which  lowers  the  emitter  sheath 
potential  drop  (|V£I,  Eq.A-3)  and  this  in  turn  lowers  the  diode  voltage  (Eq.A-l). 


Effect  of  Electrode  Spacing  (d) 

The  effect  of  the  electrode  spacing  (d)  on  J  +  and  ng  is  shown  in 


TAble  I-C  and  the  effect  on  the  I-V  characteristics  in  Fig.  k.  A  value  of 
J  >1  A  cm"2  was  obtained  for  d=0.5  cm  at  S=1.0xl0  7cm  3sec  but  again  the 


electron  density  is  beyond  the  range  oi  applicability  of  the  diffusion  model, 


Eq.(l5). 


11 


13 


VII .  REFERENCES 


C.  B.  Leffert,  D.  B.  Rees,  and  F.  E.  Jamerson,  J.  Appl.  Phys .  133  (1966) 

C.  B.  Leffert,  D.  B.  Rees,  and  F.  E.  Jamerson,  ONR  Annual  Report  No.  5, 
Contract  Nonr-3109(00)  (Oct.,  1964)  • 

C.  B.  Leffert,  D.  B.  Rees,  and  F.  E.  Gifford,  ONR  Annual  Report  No.  6, 

Contract  Nonr-3109(00)  (Oct .  1965)* 

C.  B.  Leffert,  Adv.  Energy  Conv.  2,  417  (1962). 


APPENDIX  A  -  PROBLUn  STATEMENT 

ELECTRON  TRANSPORT  PROGRAM  (DIFFUSION  ONLY) 

Internal  Diode  Voltage  from  Net  Diode  Current 

(Diode  Voltage  =  V  +  (^  -  »E) 

Solve  for  the  internal  diode  voltage  V,  where 


<  'net  <  Jnet> 


v  =  -  (v-  +  vp  -  Tp) 


and  for  (J^  <  J„et  < 


V  -  -  (Vp  -  V„  -  Vp) 


Symbols  are  defined  In  Table  A-I  ana  the  boundary  currents  Jnet,  and 


J"“T  will  be  defined  last, 
net 


The  emitter  sheath  potential  V£  is  given  by 


11,600 


ffnet  +  no  VE  | 
l  VE  +  no  VE+ 


the  collector  sheath  potential  by 
for  (^e?CJnet  ^  Jnet) 


.  .  |  TC  |  [VC  +  >net^nd  1 

\ 11,600 1  ln  „  -  J 


_  C  ^  ^  Tmax\ 

for  (,net<Jnet<  Jnet) 


11,600 


c  .  Sd-VE  lVE  +  noV 

_  In  — — r - 


fnet  +noY 


15 


”d  -  f  -  i  “here 


(A-13) 


for(iCt<  Jnet£Jnet) 


A  =  1  + 


(A-14) 


d  u+  \  P  ,  Sd2 

B  =  r  -r1-:  lnet  -  rD 

a  \i  +  [i  I  a 


(A-15) 


+  + 
+  ^  VE  VC 

0  -  V  ve  +  -r— 


(A-16) 


+  p  +  veVb  .  Vc+ 

VE  'net  “  ve  Sd  ‘  A  A 


(A-1T) 


VpV_+B 

E  .  Sd  Vg 


(A-18) 


,  C  ,  ,  _max\ 

and  for  (Jnet<  Jnet<  J„et) 


A  =  1  + 


(A-19) 


D  +  X  ' 

a  \  n  +  \x 


net  - 


Sd2 
2  D 


(A-20) 


1  +  ,  VE  VC 

VE  VE  A 


(A-21) 


D  =  - 


+  VE  VC  B  L  VEVC 
Sd  VE  '  A -  +  "A 


(A-22) 


„  Vc B 

E  =  -  J 


‘  Sd  VE  •  Hiet  VE 


(A-23) 


The  random  particle  current  density  V?  of  the  charge  cairiers 
(j  =  :  for  electrons,  =  +  for  ions)  in  the  plasma  near  the  electrodes 
(i  =  E  for  emitter,  =  C  for  collector)  are  defined  in  terms  of  the  particle 
density  (n=n”=n+)  and  the  average  velocities  v^  as 


where 


C3  - 


n^ 

“T" 


T 

8  k  i 
n  nu 


1/ 2 


(A-2U) 


(A-25) 


We  use  the  rationalized  MKS  system  except  that  we  prefer  the^number  densi¬ 
ties  in  units  of  cm-3  and  the  velocities  in  units  of  cm  sec'  and  since  the 
number  1+  always  appears  with  v^  ,  we  set 


n,cm 


x  , cm  sec 


-1 


(A-26) 


where 


-1 

cm  sec 


(A-27) 


This  equation  is  used  to  calculate  v£  ,v£  ,vc  and  vc  . 

The  mobilities  of  the  charge  carriers,  nJ,_  are  defined  in  terms  of  the 
mobilities  at  1  torr  and  273 °K,  by 


J  = 


(A-28) 


where  p  is  the  gas  pressure  in  torr  and  j  is  defined  above. 

The  minimum  current  density  (J^)  is  the  value  of  the  current  density 
as  the  collector  is  driven  negative  with  respect  to  the  emitter.  It  is 
given  by 


Tmin 

net 


min  „  min  h 
net  =  ‘  4  "d  VC 


4  v. 


nmln+  P 
o 


(A-29) 
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where 


(A-30) 


f>  =  (Sd^)/2  Dp 


and 


1  + 


D  1  +  - 

a  \|i  +  [i 


(A-31) 


4  -I  £ 

B  B 


and  >sin-(§ti[-i+V 

The  factors  in  (A-3$  are  given  by 

A  =  VE  VE  > 


B  =  (vevc+/5)  -  Sd  vE+ 


(A-32) 


(A-33) 

(A-34) 


(vePvc+/5)  -  Sd  vw 


E 


(A-35) 


The  maximum  current  (jraax)  is  the  maximum  value  of  the  net  current 

net 

density  or  the  saturation  value  of  the  current  density  as  the  collector  is 
driven  positive  with  respect  to  the  emitter.  It  is  given  by 


max  P max  =  q  nmax 


'net  ^  ’net 


'l  nd  VC  =  q  VC 


nmax  +  p| 
o 


(A-36) 


where  |3  is  given  by  (A-30)  but  here 


and 


where 


and 


- 1  +'\Ji  -(HIST. 

A  =  VE_VE+  +  (vcfvE+/&)’ 

B  =  (P  vc've+/5^  _  Sd  VE+ 

C  =  -  Sd  vE 
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max 


n 


B_ 

2A 


(A- 37) 

(A-38 

(A-39) 

(A-Ij-O) 

(A-41) 


and  p  is  given  by  (A-30) . 


TABLE  A-I .  Definition  of  symbols  ♦ 


SYMBOL 


Fixed  Constants 


FORTRAN  NAME 


Input  Constants 
d,  cm 

e  -3  -1 

S,  cm  sec 

-  2  -1 

Da, cm  sec 

n  2-1 

D+,cm  sec 

.  -2 
Vg,  A  cm 

kg 

+  2  ..-1  -1 

li0,cm  volt  sec 

-  2  ..-1  -1 

H0,cm  volt  sec 

te-'k 

TC'°K 

V 1 

p.  torr _ 

I ndependent  Variable 

J  ,  A  cirf^ 
nett _ 

Dependent  Variable 
V,  volt 

V£,  volt 

Vc ,  volt 
Vp,  volt 

Jmi"  A  cm'2 

.2 

Jnet>  ^  cm 

JneL  cm"2  sec_1 

P net >  cm  2  sec 

pmin  ~2  -1 

|  lnet,  cm  sec 

r  C  -2-1 

I  net*  cm  sec 
nmax  _2 

[  net,  cm  sec  1 

n(x),  cm‘J 
n0,  cm- 3 

na,  cm-3 
vg+ ,  cm  sect 
vg”,  cm  sec" 
vc+,  cm  sec~l 

VC  ,  cm  sec-'*' 


JRICH 


MU0I6N 

MUOEL 


JNET 


VAPPL 

VOLTE 

VOLTC 

VOLTP 

JNET(NJ+l) 

JNET(NJ+2) 

JNET(NJ+3) 

GAMMA 

GAMMIN 

GAMMAC 

GAMMAX 

RNX 

NO 

ND 

VEION 

VEEL 

VCION 

VCEL 


DEFINITION 


Electron  charge 


Boltzmann  constant 

Mass  of  electron _ _ 

Spacing 

Ion  generation  rate 
Ambipolar  diffusion  coefficient 
Ion  diffusion  coefficient 
Richardson  emitter  current 
Mass  of  ion 

Mobilities  at  1  torr  and  273°K 

Emitter  temperature 

Collector  temperature 

Mean  gas  temperature 

Gas  pressure _ 

Net  diode  current 


VALUE 


1 .6021xl0'19coul . 

3-14159 

1.3804x10-2V°k 

9.1Q84xlO~31kg 


1 .10x10 
1 .20X101 
6  .00x10° 
2.00 

6 .6801x10' 
5-32X103 

5 .ooxio6 

1700. 


1300. 

1.0  x  102 


1.0  x  10' 


Diode  voltage  (for  <t>c=$E) 

Emitter  sheath  potential  difference  Vp=((^V)p| 

Collector  sheath  potential  difference  Vq=\  ( A  V)c | 
Plasma  potential  difference  Vp=  |(A.v)P| 

Saturation  back  current  density  (V  - >-©o) 

Current  density  at  which  Vq=0 

Saturation  forward  current  density  (V  - >+00) 

r net  =  Jnet/^ 
rmin  _ 

I  net  -  Jngt/^ 

r  net  =  Jnet/ Q. 
rmax  _  Tmax/ 

I  net  “  Jnet'<l 

Plasma  density  at  a  distance  x  cm  from  the  emitter 
Plasma  density  near  emitter  sheath 
Plasma  density  near  collector  sheath 
Average  velocity  of  ions  in  plasma  near  emitter/4 
Average  velocity  of  electrons  in  plasma  near  emitter/4 
Average  velocity  of  ions  in  plasma  near  collector/4 

Average  velocity  of  electrons  in  plasma  near  cc01ectar/4 
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APPENDIX  B 


, _ TABLE  B-I.  Program  listing  and  sample  output. 

[BFTC  MN5  LIST.REF.M94.XR7 

ELECTRON  TRANSPORT  PROGRAM  (DIFFUSION  ONLY) 


JOB  NO 


calculates'  internal  diode  voltage  (VAPPL  b“^lt^"‘°^^^C+V0LTP) 

FOR  INPUT  VALUES  OF  NET  CURRENT  DENSITY  (JNET.  AMPS/CM**2) 

DIODE  OUTPUT  VOLTAGF  WOULD  BE  =  (VAPPL  +  CONTACT  POTENTIAL) 

USES  SUBROUTINE  VFORIS  FOR  ELFCTRON  RETARDING  COLLECTOR  SHEATH 
USES  SUBROUTINE  VFOR1L  FOR  ION  RETARDING  COLLECTOR  SHEAT 

1  NEUTRAL  PLASMA  (ION  DENSITY  •  ELECTRON  DENSITY  -  NX ( X )  ) 

2  COLLISION  DOMINATED  PLASMA  (ELECTRON  MEAN  FREE  PATH 

GREATER  THAN  DEBYE  LENGTH  BUT  LESS  THAN  DIODE  SPACING) 

3  DIFFUSION  CONTROLLED  PLASMA  (RECOMBINATION  LOSS  ■  NIL) 

4  UNIFORM  ION  GENERATION  RATE  (S> 

ASSUMES  5  COLLISION  LESS  SHEATHS 

6  ELECTRON  RICH  EMITTER  SHEATH  (ION  RETARDING  TO  PLASMA) 

7  MAXWELLIAN  DISTRIBUTION  OF  ELECTRON  AND  ION  VELOCITIES 

8  ELECTRON  TFMP  «  ION  TEMP  *  GAS  TEMP  ■  ELECTRODE  TEMPS 

9  PLASMA  TEMP  *  MEAN  ELECTRODE  TEMP 

PRINTS  1  VAPPL  VS  JNET 

2  JNET (MIN) . JNET ( MAX ) ♦ JNET  <  VC*0 ) 

3  NX(X)  VS  X  FOR  EACH  JNET  INPUT 


DIMENSION  T I TLE ( 24 ) « JNET ( 50 ) • VAPP ( 50 ) 

DIMENSION  RNO < 50 ) » RND ( 50 ) 

COMMON  /C0M1/D.S.DA.DP. JRICH.O.Pl .KB. ME. Ml .MUEL.MUlON. 

1 VEEL . VCEL . VE I  ON « VC  I  ON . KVEL « <V I ON ♦ NU . TE « TC « TM «N J » LD 

DOUBLE  PRECISION  P . MU  I  ON . MUFL «0 * P I « *’B *  ME .Ml  » TE . TC # TM « MUO I  ON. MUOEL 
DOUBLE  PRECISION  D.S«DA,DP.JRICH, JNET. GAMMA. VAPPL. VAPP. NU 
DOUBLF  PRECISION  VEEL . VCEL , VE I  ON . VC  I  ON . K VEL . KV I  ON 
DOUBLE  PRECISION  NO .NO 
DOUBLE  PRECISION  NOC .NDC . VAPPC 

DOUBLE  PRECISION  RNO . RND . BETA .DELTA . AA «BB . CC «NOM« NDM .GAMMAX 
DOUBLE  PRECISION  NOM I N .NDM I N . GAMM IN 
DOUBLE  PRECISION  ALFA » B 1 . B2 «B3 . B5 . A . B « C « GAMMAC 
NAMEL1ST/INPUT1/D.S.DA.DP. JRICH.Ml . MUO I  ON . MUOEL » TE • TC . TM . P 
NAMEL I ST/INPUT2/NJ.LD 

NAMEL I  ST/NAM 1 /K VEL .KV I  ON  « MU  I ON.MUEL  * NU . VEEL . VCEL . VE ION. VC  I  ON 
NAMEL I ST/NAM2/N0M.NDM. GAMMAX 
NAMEL I ST/NAM3/N0M IN. NDM IN. GAMM I N 
NAMFL I ST/NAM4/N0C « NDC  » GAMMAC .VAPPC 


PO  RFAD  <5»100)TITLE 
100  FORMAT  <  1 2A6) 

I  OB  WRITE (6.1 10)TITLE 

I I  0  FORMAT  (1H1.20X.12A6/21X.12A6) 
120  READ  ( 5. INPUT  1 > 

125  WRITE (6. INPUT  1) 

130  DO  135  L«1 .50 

131  JNET ( L  )  =0 .000 

134  RNO( L ) *0.000 

135  RNO <L) *0.000 
140  RFAD  (5.INPUT2) 

145  WRITE <6. INPUT?) 


22 


1-5(6 

150 

151 
160 
161 

170 
175 
160 
166 
190 
196 
200 
P  05 
210 
220 
225 
230 

235 

236 

240 
250 
260 
_  270 
~  260 
290 
300 
310 
320 
330 
340 
350 


NJ3=NJ+3 

RF  AO  (5*151  ) (JNET(N) *N 
FORMAT  <  60 12*6) 

WRITE  (6. 161 ) < JNET(N) *N 
FORMAT (1H  «  27H I NPUT  FOR 


1  «NJ> 

= 1 «NJ) 

NET  CURRENT* 


AMPS/ ( 5E20 .8  > ) 


Q  *  1 .60210-19 
PI «3. 1 4159  DO 
KB»1 .36040-23 
MF*° • 1 0640-31 

KVFL  =1 .0D2*DSQRT< ( KB/ME >#< 0 • 5D0/PI 5 ) 

KV10N=1 .002*OSOPT  ( (KB/MI )#<0*5D0/Pt >  > 

MU  I ON*MUO I  ON/P 

MUFL  =MUOFL/P 

NU*  JR  I CH/Q 

VEEL=KVEL*DSQRT  <  TF ) 

VCEL=KVEL*DSQRT  <  TC ) 

VE I ON=KV I ON*DSQRT ( TE ) 

VC  I ON  =  KV I ON*DSQRT ( TC ) 

WRITE (6.NAM1 1 


BFTA*S#D*D/(2.0D0#DA)  — 

OFLTA=  ( 1 .000+1 D/OAHMMUI  ON/  (MUION+MUEL)  >*VCEL> 

AArn ( VFFL*VF 1 0N+ ( VCEL*VE I  ON  1 /DELTA  1 

rr«  ( BFTA/DELT  A )  *VCEL*VE  I  ON— S*D*VE  I  ON  . . 

rr,  000#  AA  n  #  (  - 1  .  ODO+DSQRT  (  t  .  ODO-4.0DO*  (  AA/BB )  *  ( CC/BB )  >  ) 

NDM* ( NOM+BETA 1 /DELTA  -  •  . 

GAMMAX*NDM#VCEL  .  -  . 

WRITE(6«NAM2)  -  - . - . 

RN0(NJ+3> =NOM  -  '• 

RND(NJ+3) =NDM 
JNET < NJ+3 ) =0*GAMMAX 


360 

362 

364 

366 

368 

370 

372 

374 

376 

37P 

3B0 

400 
4  1  0 
ApO 
430 
440 
460 
460 
470 
460 

462 

463 
490 
600 


DELTA  =< 1 . ODO  +  <  D/D A 1 #  <  MUEL/ « MU I ON+MUEL  > )#VCION) 
aa=vffl*veion 

BB* ( NU/OFLT  A ) *VC I  ON— S#D#V£ I  ON 

CC*NU* ( BET  A /DELTA ) #VC 10N-S#0#NU  _  #  a  a  saa \ m  t CC/BB 111 

N0MIN= (BB/( 2.0D0#AA) > # < - 1 • ODO+DSQRT < 1 .000-4# 000* ( AA/SB)*( CC/BB) 1 

NDMIN= (NOMlN+BETAl/DELTA 
GAMM I N  =  -NDM I N*VC I  ON 
WR I TE ( 6  «NAM3 1 

RN0(NJ+1 > =N0MIN 

RN0(NJ+1 ) =NDM1N 
JNFT ( NJ+1 ) =0*GAMMIN 

Al  F A  = ( 1 . ODO+VC 1 ON#D/D A  > / ( VCEL- VC  I  ON  M-  <  MU  I  ON/ 1  MU  I ON+MUEL >  >  *  < D/DA ) 

B 1  =VFFL 
B2«VF 1  ON/MU 
B3=S*D 

B6«VCFL/< VCEL-VCION) 

B  I- 1  !ooO-k?«3l:ALFr‘B5*4LF-2  *000.61  .B2.ALFA.0ETAtB5-B2.B5.DETA 

C  «•  R2#03#BETA  +Bl*B2*BETA*BETA 

GAMMAC=(R/<2.0D0*A) )*(-!. ODO-DSORT ( 1 . 0D0-4. ODO* <  A/B 1  *  <  C/B  >  >  > 

FORMAT( 1H  ?24HCALCULAT10N6  FOP  GAMM AC*  *E16.8« 15H  WHERE  VC  IS  0) 
CALL  VFOR I S  <  GAMM AC  *  NO ♦ ND ♦ VAPPL ) 

NOC=NO 
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*510  NDC=ND 

520  vappc=vappl 

530  WRITE  (6.NAM4) 

540  RNO(NJ+2> =N0C 
550  RND(NJ+2> =  NDC 
560  VAPP(NJ+2 ) =VAPPC 
570  JNET(NJ+2)=0*GAMMAC 

600  DO  BOO  Nl=l .NJ 

610  GAMMA=JNET(N1 I/O 

614  IF  <  GAMMA-GAMM IN ) 616.616.620 

616  VAPP ( N 1 )=-l *0010 

618  GO  TO  000 

620  I F ( GAMMA-GAMMAX ) 650  *630 »  630 
630  VAPP ( N 1 >=1 .0010 
640  GO  TO  800 

650  I F  <  GAMMA-GAMMAC ) 660 ♦ 660 »  680 

660  WRITE  (6.661)  JNET(Nl) 

661  FORMAT ( 1H  « 23HC ALCULAT I ONS  FOR  JNET  «.E16.B.57H  WHERE  VC  IS  -VOLTC 
6621 (IE.  ELECTRON  RETARDING  TO  THE  PLASMA)) 

665  CALL  VFOR I S  < GAMMA • NO t  ND • VAPPL ) 

670  GO  TO  700 

680  WRITE  (6.681)  JNET  ( Nl  ) 

681  FORMAT ( 1 H  « 23HC ALCULAT I ONS  FOR  JNET  =»E16.8.57H  WHERE  VC  IS  +VOLTC 
6821 (IE.  ION  RETARDING  TO  THE  PLASMA)) 

685  CALL  VFOR I L < GAMMA « NO »ND. VAPPL ) 

700  RNO ( N 1 ) *N0  , 

710  RND  <  N 1 )=ND 
720  VAPP(N1 ) » VAPPL 
800  CONTINUE 

900  WRITE  (6.910) 

910  FORMAT  ( 1H1 « 13HFINAL  SUMMARY  /1H0» 13X.4HJNET  *26X.2HN0«25X.2HND« 

91 1 124X.5HVAPPL) 

920  WRITE  (6.930 ) ( JNET ( I ) «RNO(  I  ) «RND( I > .VAPP  <  I  )  *  I  ■  1 *NJ3) 

930  FORMAT  (4028.16) 

950  GO  TO  20 
END 

. IBFTC  SUR6  LIST.REF.M94.XR7 

SUBROUT  I NE  VFOR I S ( GAMMA .  NO  ♦ NO . VAPPL ) 

COMPUTATES  INTERNAL  VOLTAGE (VAPPL)  FOR  A  GIVEN  NET  CURRENT  (GAMMA) 
FOR  COLLECTOR  SHEATH  ELECTRON  RETARDING  TO  THE  PLASMA  (VC*-VOLTC) 

ALSO  RETURNS  PLASMA  DENSITY  AT  EMITTER  (NO)  AND  COLLECTOR  (NO) 

FOP  EACH  CALL  THE  PLASMA  DENSITY  (Rmx)  IS  PRINTED  OUT  FOR  LO  STEPS 
OF  DISTANCE  (X)  ACROSS  THE  GAP  (D) 

ALSO  PRINTED  ARE  VOLTE.  VOLTC*  AND  PLASMA  DROP  (VOLTP) 


COMMON  /COM  1 /D .S.OA.DP. JR  I CH«  O.PI .KB.ME.Ml *MUEL*MUlON» 

1 VEEL . VCEL . VE I  ON . VC  1  ON « KVEL  »KV I  ON. NU «  TE . TC . TM .NJ.LD 
DOUBLE  PRECISION  P » MU  I  ON  .  MUEL  .O.PI  .  ►'B  .ME  .Ml  .TE.TC.TM 
DOUBLE  PRECISION  D « S . DA . OP « JR ICH. JNET . GAMMA . VAPPL .VAPP. NU 
DOUBLE  PRECISION  VEEL . VCEL . VE I  ON. VC  I  ON. KVEL* KV ION 

DOUBLE  PRECISION  NO.ND 

DOUBLE  PRECISION  X.RNX 

DOUBLF  PRECISION  AN.BN.CN.DN.EN 


2k 


OOUBLE  PRECISION  AVP . BVP . CVP . KVE , KVC , VOLTE .VOLTC . DVP l . DVP2 
DOUBLE  PRECISION  VP1 « VP21 .VP22.VP2.V0LTP.VAPPL 
NAMEL I  ST /N AM6/N0 . ND 
NAMEL I ST/NAM7/KVE  « KVC  *  VOLTE  *  VOLTC 
NAMEL I ST/NAM8/V0LTP . VAPPL 
C 

100  AN»1«0D0  +  D*VC I  ON/DA 

1  1  0  BN* ( ( D*MU I ON*GAMVA ) / ( DA*  <  MU  I ON+MUEL  >  >  ) -S*D*D/ ( 2  »ODO*DA ) 

120  CN*VEEL*VE I ON  +  VE I  ON* VC  I  ON/AN 

1 30  DN*VE I ON*G AMMA-VE I ON*S*D- VE I ON*VC I ON*BN/AN+NU*VC 1 ON/AN 
140  FN*-( ( (NU/AN)*VCION*BN)+(S+D*NU> 1 

160  NO*— ( ( 0*5DO#DN/CN)*  < 1 . ODO-DSORT ( 1 »ODO-(4.0DO*(EN/DN) /(DN/CN) ) ) ) > 

170  ND*NO/AN  -  BN/AN 

180  WR I TE  (  6 « N AM6 1  ..  _ _  - 

200  KVE*  (  GAMMA+NO*VEEL  )  /  (  NU+NO*VE  I  ON)  .  . . . 

210  KVC*  (  VC  1  ON-fGAMMA/NO  ) /VCEL  . . . . . - .  . 

220  VOLTE  =  -(  TE/1  •  16D4)*DL0G<  KVE)  _  _  _ _ _ 

230  VOLTC*  —  <  TC/1  •  1604 )  *DLOG(  KVC  )  .  .  . . 

240  WRITE  (6.NAM7)  _  . . . 

250  XsO.ODO  ,  . .  . 

252  XD*LD 

254  XDD*D/XD  . J . .  . 

255  LD1*LD+1  ...  . . - . - 

256  WRITE  (6.257)  GAMMA 

257  FORMAT (  1 H  « 18HNX  VS  X  FOR  GAMMA* ,E 16.8//  1 3X* I HX. 27X« 2HNX) 

260  DO  290  L* 1 *LD 1 

270  RNX*NO—  (  S/(  2.0D0*DA  )  )*X*X-(  (MUlON/lMUlON-fMUEL)  )* (GAMMA/DA)- . 

""  271 1NO*VEION*KVE/DA)*X 

. 275  X*X+XDD  .  . . . . 

280  WRITE  (6.281)  X.RNX  .  .  „  .  . 

281  FORMAT (2D28. 16)  , 

290  CONTINUE 

300  AVP*  <  S*D— ND*VC I  ON-MU  I ON*GAMMA/ (MU I ON+MUEL ) ) /DA 
310  BVP*2.0DC*AVP*DA/S 
320  CVP*2.0D0*N0*DA/S 

400  DVP 1  * ( DSORT ( BVP*BVP+4 . 0D0*CVP ) -BVP ) /2 • ODO 
410  DVP2* (DSORT (BVP*BVP+4.0D0*CVP)+BVP)/2« ODO 
420  VP 1  * ( DP/MU  I  ON ) * ( DLOG  <  C VP+BVP*D-D*D  > -DLOG ( CVP ) ) 

430  VP2 1  * ( 4 • ODO*DP*GAMMA ) / ( S# ( MU I ON+MUEL  >  *DSQRT ( BVP*BVP+4  *ODO*CVP ) ) 

440  VP22*DLOG( (D+DVP1 )/(DVP2-D) )-DL0G(DVPi/DVP2) 

500  VP2=VP21*VP22 
510  V0LTP*VP1+VP2 
520  VAPPL=-( VOLTE+VOLTC-VOLTP) 

530  WRITE  (6.NAM8)  . .  .  . 

600  RETURN  . 

END 

%  I BFTC  SUB7  LIST.REF.M94.XR7 

SUBROUT I NE  VFOR I L  < GAMMA ♦ NO  «  ND ♦ V APPL ) 

C 

C 

C  COMPUTATES  INTERNAL  VOLTAGE ( VAPPL >  FOR  A  GIVEN  NET  CURRENT  (GAMMA) 

C  FOR  COLLECTOR  SHEATH  ION  RETARDING  TO  THE  PLASMA  <VC*+VOLTC) 

C 

C  ALSO  RETURNS  PLASMA  DENSITY  AT  EMITTER  (NO)  AND  COLLECTOR  (ND) 

C  FOR  EACH  CALL  THE  PLASMA  DENSITY  (RNX )  IS  PRINTED  OUT  FOR  LD  STEPS 

C  OF  DISTANCE  (X)  ACROSS  THE  GAP  CD) 

C  ALSO  PRINTED  APE  VOLTE.  VOLTC.  AnD  PLASMA  DROP  (VOLTP) 

C 

C 
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COMMON  /COMl/D.S.OA.DP.JnICH.O.Pt .KR.ME .Ml .MLEL.MUION. 

I VEEL .  VCEL  .  VE  ION .  VC  I  ON .  KVEL  4KV I  ON  4  NU .  TE  4  TO  4 TM  4N  J  t  0 

nmiRiF  PRFCISION  P  ,  MU  I  ON  « MUEL  , 0  ,  P  I  ,  KB  «  ME  *  M  *  *  TE  *  TC  *  „ 

DOUBLE  PRECISION  ^;^iDp4jR|CHtJNFTfGAMMA*VAPPL.VAPP,NU 

VEEL  , VCEL ♦ VE I  ON  «  VC  I  ON « K VEL ♦ KV I ON 
NO  «ND 
X.RNX 

AN.BN.CN.ON.EN 

AVP.BVP.CVP.KVE  ,  KVC, VOLTE. VOLTC .DVPl *DVP2 
VPl  ♦ VP21 « VP22 « VP2  *  VOLTP « VAPPL 


OOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

double 


PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 


•  NAMFl  l ST/N AM6/NO  t  NO 

NAMEL I ST/N AM7/KVE ♦ KVC « VOLTE « VOLTC 
NAMEL I ST/NAMB/VOLTP ♦ VAPPL 

110  BN«-nD*MUEL#GAMMA)/<DA*<MUlON+MUELn  >-S*D*D/(2»ODO*DA > 

120  CN=VEEL*VEION+VEION*VCEL/AN 

1 30  DN*-VE I ON*S*0-VE I ON*VCEL*BN/AN+NU*VCEL/AN 

i:°0  no:;'<  !^^:^cLN;:u?»o^:opTu:r.«*.oDo..EN,oN,,«DN,cN,  „ . , 

170  nd*no/an  -  bn/an  —  -  -  • 

180  WRITE ( 6.NAM6) 

200  KVE* ( GAMMA+NO*VEEL > /<NU+NO*VE I  ON) 

210  KVC«<S*D-N0*VEl0N#KVE)/(ND*VC10N)  . 

VOLTE*- (TE/1* 1604) *DLOG( KVE) 

VOLTC*- ( TC/ 1 • 1604) *DLOG( KVC ) 

WRITE (6.NAM7) . 

X“0»0D0  ...... 

XD*LD 
XDD*0/XD 
L01 *L0+1 

WRITE  (6,257)  GAMMA  ^  GAMMA*, El 6* 8//  13X. 1HX.2TX.2HNX) 


220 

230 

240 

250 

252 

254 

255 

256 

257 
260 
270 


(  (MUEL/CMUION+MUELI )*( GAMMA/DA )- 


280 

281 

290 

300 


FORMAT (1H  « 1 flHNX  VS  X 
DO  290  L* 1 »LD1 

c,..  RNXrNO-(S/(2.0D0*DA))*X#X4 

2711 (NU*(KVE/DA)-NO*<VEEL/OA) ) )*X 

275  X*X+XDD 

WRITE  (6,281)  X.RNX 
FORMAT (2D28. 16) 

A VP* ( ^NO* VE I ON*KVE  )  /DA- ( MU  I ON*GAMMA ) / <  DA* (MU I ONTMUEL  > ) ) 

310  BVP*2.0D0*AVP*0A/S 
320  CVP*2.0D0*N0#DA/S 

400  DVP1 * ( DSORT ( BVP*BVP+4 . 0D0*CVP ) -BVP ) /2 • ODO 

410  DVP2*(DSORT(BVP#BVP+4.0DO*CVP)TBVP)/2*ODO 

a»a  \/d i  /  opyML j  1  ON )  ^  ( 0l»06  ( CVP+BVP^O^D^D ) -*Dl*OG  ( CVP )  ) 

430  VP21 * ( 4 . ODO#DP*GAMMA ) / <  S*  <  MU I ON+MUEL ) *DSORT  <  BVP*BVP74 . ODO#CVP ) ) 
440  VP22*DL0G ( ( D+DVP1 > / (DVP2-D) )-DL0G(DVPl/DVP2) 

500  VP2*VP2 1 * VP22 
V0LTP=VP1+VP2 

VAPPL*- (VOLTE-VOLTC-VOLTP)  .  - 

WRITE (6.NAM8) 

RETURN  . - 

END 


510 

520 

530 

600 
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I 


TABLE  B-EE  .  Sample  of  code  output . 

RUN  10.1  V  VS  I  NE-AP 

JNET*0*03  =  1  .0,2.5,5.0,10.0,70.0,120.0 


r 


I INPIT l 


.  0 

a 

0.99999944999999970- 

01. 

E  ic .  ” 

m 

O.llGOCCOCOOOCOOOOD 

17, 

3JT 

a 

0. 12CC000CCC  CCOOCOD 

02, 

1  “OP 

a 

0.555 99 55559 99 9998D 

01, 

17  JR  ICE 

a 

C. 2C  CCCCCCCCOCOOOOD 

01, 

Em#  | 

a 

C. 668005555 9 999997D- 

25, 

l  MU0ICN 

a 

0.53 159595999999 99 D 

04, 

Mj, 

PUCEl 

m 

C.5CCCOOCCCCCCOOOOO 

07, 

0  te 

a 

0. 1  7C0GCCC00CCGCC0D 

04, 

TC 

a 

0.  8555955955999998D 

03, 

IS  ™ 

a 

0.13CCCGCCCC 0000000 

04, 

P 

r 

a 

c.icccocccccocooooo 

03, 

L  $  ENC 

•*  S1NPLT2 

L  -  - 

NJ 

a 

6 » 

[ 

a 

10, 

t  ENC 

r  INPIT 

EUR  NET  CURRENT .  AMPS 

1  0 

.45595959E- 

-03  C.25CCOOOOF- 

-02 

•  0 

.999999 99 E 

-02  0.69999999E- 

■01 

T  SNAHl 

L 

KVEL 

a 

0. 1563C72312065078D 

06, 

i»  KVICN 

* 

C. 57348391* 65  30672D 

03, 

PUICN 

fm 

a 

0.53155559595999980 

02, 

l  MUEL 

* 

C.  5CCCCCCCCC  CCOOCOD 

05, 

m  NU 

* 

C.1246361525497784D 

20, 

•  VEEL 

a 

C.64034511R686654fiD 

C7, 

|  .  VCEL 

a 

0.46552169351952340 

07, 

VEICN 

3 

0.23fc45347625535llD 

05, 

I 
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0. 499999 99E-02 
C.12C0CQ00E  CO 


\V  -  ■•  " 


VCICN  *  0. 172C451749559202D  05, 

$  ENC 
SNAH2 
NOM 
NOM 

GAMMAX 
S  ENC 


0 .2 C 8 ?54  165034 73 7 ID  12 T"'  _ _  _ 

0.33228222  146868610  1 1  *  _ _  . 

0.57167552!  54  06 1 3  20D  15, 

C ALCULAT I CNS  FOfc  GAMMAC*  C.12028628E  18  WHERE  VC  IS  0 


SNAP* 

-•  - 

NO  . _> 

C. 2 23C88 40014703000  12, 

NC  ...  *  .... 

0.25912527802746860  11, 

$  ENC 

. 

0. 12401651525831 16D  CO, 
0. 55599955999999570  00, 
C.3C59C33542201612U  00, 
0.3187105751725G70D-15, 

$  ENC 

NX  VS  X  FOR  GAMMA*  C.12028628E  18 


0.223G8840C14703000  1? 
0.61587080919670410  12 
C. 91698655354325000  12 
0.11264356331866680  13 
C. 12442 18048126957D  13 
0. 12703337583641 18D  13 
C. 12047828838981520  13 
0. 10475653C47290560  13 
0. 798681C6085683270  12 
0.45813015228148110  12 
0. 259 12 579CC3C 00970  11 


INAM8 


X 

C.5SS59998528579C30-C2 
C. 155559557E5755810-C1 
C. 29555995678653 7 10— C 1 
C. 39955995 57 159 16 10-C1 
C.  455 5 599546  4485520- Cl 
C. 5 555999935 736 7410-01 
C. 65 55 59 95 2 5 02 8 5 310- Cl 
0.75599999143183220-01 
C.65555995C36081110-C1 
C. 555559985? e975C30-Cl 
C. 1C55559988218765D  00 


$  NAM  7 
KVE 

K  VC  * 
VOLTE  * 
V0L1C  * 


INAN3 
N0M1N  * 
FOMIN  * 
GAMP  IN  * 


0.26239276537246240  12, 
0.1  1464340304141  10D  12, 
0.53414848507359820  18, 
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VGLTP  *  C.  1  7C338374  14457550  00, 

VAPPl  =  -0.  ]  -55fc45  8007  5081  50  00, 

S  ENC 
SNAP4 

NCC  *  0.22306840C1470300D  12, 

NCC  *  0.2*912527802746860  11. 

GAMMAC  •*  C.12028627586C1466D  18, 


VAPPC  *  -0.13556498CC165815D  00, 

9 


%  ENC 

CAL  COL AT ICNS  FOR  JNFT  = 

SNAP  6 


0. 55599999L-03  WHERE  VC  IS  -V0LTCI  IE, 
ELECTRON  RETARDING  TO  THE  PLASMA) 


NC  *  0.2096586393 92 39 2CD  12, 


NO  *  0.228161989C653208D  11, 


S  EM 
SNAP  7 

KVE  *  0. IC80C16902 1321380  CO, 
KVC  *  0.44515977032  79608D-01, 
VOLTE  *  0.3 26 166 7485 481 5921)  00, 
VOL TC  *  0.2414^,106963985140  00, 


$  ENC 

NX  VS  X  FUR  o ANM A  =  0.6241  EC76E  16 

X  N* 

0.599959989265/5030-02  0. 20965H6393923920D  12 

0.19595999765755810-01  0.60447439160613090  12 

C. 2555599S67869271D-CI  C. 9076? 34791 1674 180  12 

0.35555995571591610-01  0.11191059019242240  13 

0.455559994(448552!)— Cl  0.  12  389216600285790  13 

C. 5555599935738 7 mIO- 01  0 .  1 2 67C70 7 5 34 298050  13 

C. 65  555955250265310- Cl  C . 1 303553 1 82 1279030  1  3 

C.  75595999  1431  83220-01  0.13483685461228  720  13 

C. 85 99 5905c 36061 110-01  0.831*1804541471350  V? 

C.5S555598S26975C30-C1  C . 4 5 300048300 342 64 0  12 

C. 105959556821 87650  CO  0. 32816249 3890 11 7?n  11 


SNAPE 

VOLTE  -  -0.1(814872662790620  00, 

VAPPL  *  -0.75575654501591670  00, 


S  ENC 
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.w^iiwiMgWW'iW  UtMWII  !IWIM  »WOTW IWHWH  w  Mmi  bipww 

1 


jjOJE^.  Output  for  intermediate  current  values  were  deleted . 


CALCULA  T 1 CNS  FfJK  JNET  =  0 . 65999999E-01  WHE9E  VC  IS  fVHLTCIIE, 

KIN  RETARDING  TO  THE  PLASMA) 

$NAP6 

KC  =  0.25441267E4&456970  12, 

NO  «  C.53 7 993267530191 60  11, 


$  ENC 

SNAM7 

KVE 

KVC 

VCLTE 

VCLTC 


C. 165621488CP729150  OC, 
0.6 99 65277821016880- 01, 
0. 2fc36E444838C07760  CO, 
0.2 120 84 45528 99 30CO  CO, 


$  ENC 

NX  VS  X  FOR  GAMMA=  C.42652653E  18 


X 

0.5  5  5959985  2  8975  00-02 
C. 15559999 78 5755 810- Cl 
0.25555995 67 8653 7 10- Cl 
C. 35559999 57 1591 6 IU-C1 
0. 4555599946 44«5 520- Cl 
C. 55555999  35  73E74 IC-C1 
C. 655959 99 25028531 0-01 
C.  75595999 15 318 3220-01 
C.8555599SC360E1 110-01 
C.  5559599692 8975C 30-C l 
C.  1C59599568216 7650  CO 


NX 

0.2544126784045657D  12 
C. 65085133955625950  12 
0. 95 56 2333592482 130  1? 
C.  1 1687266675902  551)  13 
0.12901673345525600  13 
0.13199393)68117380  13 
0.  1 2  53044  6  74 3677860  13 
0.  11045833472207  070  13 
0.85925535537049920  12 
C.52236C 65 881716 360  12 
0.9375937756C699690  11 


$NAH6 

VCLTP  *  0. 1 1 1561 39655969360  01, 

V*ppl  =  0. 1C6A013956506838D  01, 


)  ENC 
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FIN  U  SLHM/mv _ 

JNET 

0, 1C CCOOOCOOCOOOCOD-02 
G.25CCCOOOOOCOOOCOD-02 
”  C . 50 OCCOCCCOCCOO CCD-02 
0.  lOCOOOOCCOCOOOOOO-Oi 
'C.65C.5559995999998D-01 
“  C.12COCOOOOCCOOCCOD  00 
-C. 9158813(566156390-04 
0.19  271064255554C8D-01 
""  'O. 85575928 79364115D-C1 


NO 

0. 2096584 39  392  39200  12 
0.21CR 3402334949 530  12 
0.21276171272109120  12 
0.21650562189422330  12 
0.25441267348456970  12 
0. OOCOOCCCCOOOOOOOD-38 
0.20879416503473710  12 
0.22306840C14703000  12 
0.26239216537246240  12 


NO 

0.32 81619890 65 32080  1 1 . 
0.32249932591955200  11 
0.3  1  30  5  938  710  0  0  4  8  3  0  11 
0. 2941 71 78840998860  11 
0.93799326753019160  11 
0.00000000000000000-38 
0.33228222146868610  11 
■'0.25912527802746660  11 
0. 1 1464340 30414110012 


VAPPl 

-0.75575654501591670  00 
-0.65650335616197850  00 
-0.550361  62062915900  00 
-0. 389 209 1741 3810 6R0  00 
0.10640139965068380  01 
0.10000000000000000  11 
-0. 00000000000000000-19 
-0.13556498007558150  00 
-0. 00000000000000000-19 
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